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ABSTRACT 


In  scenarios  such  as  software  agents  sharing  time  slots  or  blocks  in  differ¬ 
ent  physical  CPUs  or  memories,  communications  agents  accessing  channels  or 
bandwidth  on  shared  links  or  bearers  and  economic  agents  trading  in  com¬ 
modities  in  financial  markets,  the  agents  in  question  can  exchange  resources 
among  themselves  for  mutual  benefit.  This  report  develops  a  combinatorial 
probability  model  to  describe  the  like-for-like  exchange  of  resources  between 
multiple  unbiased  agents.  The  expected  exchange  rates  are  computed  for  in¬ 
dividual  agents,  syndicates  of  agents  and  the  collection  of  all  agents.  These 
performance  benchmarks  provide  intuition  and  understanding  for  the  perfor¬ 
mance  of  real  crossing  networks.  A  number  of  combinatorial  identities  are  also 
produced  as  a  consequence  of  the  modelling  and  analysis. 
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Combinatorics  of  Crossing  Networks 


EXECUTIVE  SUMMARY 

In  Defence  and  commercial  networks  there  is  a  trend  toward  user  agents  sharing  non¬ 
exclusive  access  to  common  distributed  resources  (e.g.  servers,  databases,  links  etc).  The 
modelling  of  the  exchange  of  such  resources  between  agents  is  therefore  a  problem  of 
interest. 

This  report  develops  a  new  combinatorial  probability  model  for  the  exchange  of  re¬ 
sources  between  multiple  unbiased  agents  (or  users)  within  an  arena  which  facilitates  such 
exchange  (a  crossing  network).  Once  formulated,  the  model  allows  the  calculation  of  the 
expected  exchange  rate  (the  mean  number  of  exchanged  resources)  when  an  arbitrary 
number  of  agents  meet.  For  only  a  few  (two  or  three)  agents,  many  more  probabilistic 
properties  of  the  exchange  rate  are  derived. 

The  model  uses  the  discrete  probability  theory  of  W.  E.  Deming  (of  post-WWII  TQM 
fame  in  Japan)  and  G.  J.  Glasser  from  the  1950s  to  build  a  model  which  both  further  de¬ 
velops  the  theory  and  gives  relevant  insights  into  modern  electronic  exchange  applications. 

The  expected  exchange  rate  plots  which  result  provide  benchmark  performance  curves, 
indexed  by  the  parameters  of  the  model,  for  applications  such  as  concurrent  processes 
sharing  memory  or  CPU  resources,  communications  agents  sharing  channel  or  band-width 
resources  on  common  links  or  bearers  and  economic  agents  with  a  shared  interest  in  trading 
financial  commodities. 

Being  a  mathematical  treatment,  the  report  is  necessarily,  compact,  brief  and  con¬ 
cise.  Although  motivated  by  physical  applications  (as  mentioned  above),  these  are  better 
discussed  in  separate  forums.  The  implications  for  these  applications  are,  however,  real, 
tangible  and  of  significance  to  users  of  the  corresponding  applied  systems. 
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1  Introduction 

A  crossing  network  is  the  generic  name  for  a  venue  or  mechanism  (real  or  virtual)  where 
the  current  ‘order  and  disposal’  states  of  agents  are  compared  and  resources  exchanged  - 
resources  are  ‘crossed’  between  agents  as  compatibilities  in  required  and  released  resources 
are  found.  Such  a  situation  is  typical  of  exchange  between  software  agents  or  exchange  in 
physical  markets. 

A  general  statement  of  the  exchange  problem  is  as  follows.  Consider  L  >  2  agents 
ai,...  ,aL  who  use  resources  only  from  their  own  respective  pools  Pi, . . .  ,  PL.  Each  pool 
Pi  consists  of  a  collection  i?i, . . .  ,  'tfjV;  of  resource  names.  Each  of  these  resources  is  divisible 
in  the  sense  that  multiple  agents  ai,...  ,at  may  simultaneously  use  a  different  piece  of 
a  resource  d{  if  its  name  appears  in  the  pools  Pu...  ,Pt.  With  each  agent  ah  there  is 
associated  a  current  state  list:  (Mi, . . .  ,  snidn,),  s*  6  {  +  1,-1},  of  length  n,  with  every 
di  belonging  to  Pt  and  occurring  at  most  once.  If  Si  =  +1,  agent  a,  requires  resource  di, 
while  if  Si  =  -1,  agent  at  releases  resource  i9j.  The  state  list  is  therefore  an  ‘order  and 
disposal’  list  representing  the  current  state  of  required  and  surplus  resources  for  agent  a,,. 
Such  a  list  is  used  by  any  resource  controller/allocator  serving  the  agents.  A  resource 
name  'di  is  exchangeable  between  agents  and  am  if  di  is  in  both  state  lists  and  in  a,’s 
list  is  of  opposing  sign  to  in  am's  list,  i.e.  di  is  simultaneously  surplus  to  and  required  by 
different  agents.  Any  system  of  agents  would  seek  maximum  exchangeability  to  maximise 
utilisation  of  the  resources  available  through  exchange  between  agents  and  thereby  avoid 
having  to  either  create  unavailable  resources  or  waste  unrequired  resources  on  offer  in  the 
state  lists. 

Consider  the  simultaneous  comparison  between  all  agents’  state  lists.  The  exchange 
problem  (or  crossing  problem )  is  one  of  performance  evaluation  which  asks  for  the  exchange 
rates  given  by: 

(a)  the  expected  number  of  resource  names  exchanged  between  all  L  agents; 

(b)  the  expected  number  of  resource  names  exchanged  within  the  state  list  of  a  specific 
agent  a f,  and,  more  generally, 

(c)  the  expected  number  of  resource  names  exchanged  within  the  state  lists  of  the  specific 
agents  1 , . . .  ,aes,  2  <  s  <  L  —  1 . 

Determining  these  rates  requires  the  calculation  of  probabilistic  expectations  as  per¬ 
formance  benchmarks  of  crossing  networks.  Clearly  the  result  is  greatly  influenced  by  the 
prior  distribution  assigned  to  the  occurrence  of  resource  names  in  state  lists.  The  analysis 
here  makes  the  unbiased  assumption  of  equal  likelihood  of  all  resource  names  and  equal 
likelihood  of  the  require/release  intention  for  each  resource.  This  ‘random’  collection  of 
agents  might  be  called  typical  of  non-cooperative  (independent)  agents  and  so  provides  an 
unbiased  benchmark  for  the  crossing  network’s  performance. 

This  report  gives  an  account  of  the  mathematical  detail  of  the  combinatorial  probability 
model  of  the  exchange  of  like-for-like  resources  among  many  unbiased  agents.  While 
the  crossing  problem  considered  here  is  motivated  by  a  concrete  real-world  application, 
the  formulated  mathematical  and  probability  model,  although  an  abstraction,  provides 
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concrete  exact  results  from  which  intuition  on  the  behaviour  of  real  crossing  networks  can 
be  gained. 

Under  the  unbiased  assumptions,  some  analytical  results  are  available  for  general  L. 
for  L  =  2  and  3  the  means  in  (a)  and  (b)  are  computed  along  with  some  higher  moments. 
For  arbitrary  L,  only  the  means  are  given  and  it  remains  open  to  characterise  the  higher 
moments  of  the  random  variable  describing  the  number  of  exchangeable  resource  names. 
Even  so,  the  derivation  of  the  means  and  some  higher  moments  produce  some  interesting 
combinatorial  identities  along  the  way. 

Some  specific  applications  of  the  crossing  network  scenario  are  as  follows: 

•  Each  agent  is  a  concurrent  software  process  which  may  request/release  time  slots 
on  different  CPUs,  memory  blocks  of  different  physical  memories,  or  record  blocks  in 
different  databases  (•#,); 

•  Each  agent  a;  on  a  communications  network  may  require/release  dedicated  channels 
or  capacity  on  different  physical  or  logical  links  or  bearers 

•  In  market  exchange,  if  price  discovery  is  efficient,  the  agents  a;  may  exchange  (trade) 
like  commodities  ($2)  among  themselves. 

1.1  Notation  and  definitions 

Let  ■■■it,  2  <  t  <  L.  be  the  multi-index  designator  for  the  (distinct)  agents 
ah  5  •  •  •  >au  •  A  multi-index,  being  a  combination  from  a  subset  of  agents,  never  has  any  two 
entries  equal.  Summations  with  multi-index  subscripts  are  therefore  taken  over  all  combi¬ 
nations  from  a  specified  agent  subset.  When  considering  such  a  subset  Is  =  {ii, . . .  ,is}, 
write  a/s  for  the  corresponding  agents  and  subscript  associated  quantities  with  Is.  If  s  =  1, 
denote  Is  =  {£}  by  simply  £ ,  while  if  s  —  L  (all  agents)  we  drop  the  subscript  entirely.  For 
Is  fixed,  the  agents  aI$  are  called  a  syndicate  when  2  <  s  <  L. 

Let  Ch...it  be  the  number  of  resource  names  common  to  each  of  the  pools  Pq , . . .  ,  . 

To  assist  in  visualising  results,  it  is  sometimes  useful  to  collapse  these  many  parameters 
through  the  following  functional  dependence. 

Definition.  The  L  agents  a\ , . . .  ,  are  said  to  have  exponential  decay  of  common  in¬ 
terest  (edci)  if  Ci =  ptC,  for  all  2  <  t  <  L  with  C  and  0  <  p  <  1  fixed  constants,  for 
every  multi-index  i\  •  •  •  it. 

If  0  <  p  <  1,  EDCI  conveniently  diminishes  the  size  of  pool  overlaps  with  an  increasing 
number  of  agents  considered. 

Let  A [  :=  ni/Ni  be  the  fraction  of  pool  Pi  represented  in  the  state  list  of  agent  a/. 
A  collection  of  L  agents  is  described  by  2  (f)  overlap  parameters  (Cq.-iJ,  L  list  sizes 
( ni )  and  L  pool  sizes  ( N[ )  for  a  total  of  2L  +  L  —  1  parameters.  Define  also  the  quantity 
Kk(n,  ■  ■  ■  ,it)  ■=  n^iK  —  k)/(Nir  —  k)  so  that,  for  example,  7r0(ii, . . .  ,4)  =  A q  •  •  •  \{k. 

Definition.  The  resource  name  di  is  said  to  be  uniquely  common  to  the  state  lists  of 
agents  aq , . . .  ,ait  if  it  appears  in  the  state  lists  of  these  agents  and  no  others. 
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Denote  by  the  number  of  names  common  to  the  state  lists  of  agents  a*,,. . .  ,al( 

and  by  div..ik  the  number  of  names  uniquely  common  to  the  state  lists  of  agents  ail ,ait. 
The  utility  of  the  concept  of  unique  commonality  is  that,  since  the  partition  of  common 
names  among  state  lists  into  the  uniquely  common  names  among  state  lists  is  a  disjoint  one, 
cumulative  counting  of  uniquely  common  names  avoids  double  counting  when  determining 
the  number  of  names  common  to  two  or  more  agents  within  a  specified  subset  of  all  agents. 

Let  J  (respectively,  J^,  Jjs)  be  the  random  variable  of  the  number  of  exchangeable 
names  within  the  state  lists  of  all  L  agents  (respectively,  agent  ap,  the  syndicate  ajs )  and 
let  the  mean  be  p  =  E[J]  (respectively,  fit  =  E[Je ],  p/s  =  E[J/ J).  Let  x  (respectively,  xe, 
Xfs)  be  the  random  variable  of  the  number  of  names  within  the  state  lists  of  all  L  agents 
(respectively,  agent  at,  the  syndicate  ajs )  which  are  common  to  two  or  more  state  lists 
and  let  its  mean  be  x  =  E[x]  (respectively,  Xi  =  E[xg\,  xis  —  E[xj J). 

Definition.  The  efficiency  of  the  crossing  network  (respectively,  agent  a^,  syndicate  a/s) 
is  e  :=  n/x  (respectively,  et  :=  m/xt,  e/,  :=  HiJxi,)- 


The  efficiency,  which  measures  the  ratio  of  the  expected  number  of  exchangeable  names 
to  the  expected  number  of  common  names,  gives  expected  performance  as  a  fraction  of 
the  best  expected  performance  possible. 

Let  [c]fc  :=  c(c—  1)  •  •  •  (c—k  + 1)  be  the  descending  factorial  and  (c)k  :=  c(c+l)  •  ■  •  (c  + 
fc  —  1)  the  ascending  factorial.  The  factorial  moments  of  a  random  variable  V  are  denoted 
by  ffik]  -^[[KU]-  bet  [;]  and  {^}  denote,  respectively,  the  Stirling  numbers  of  the  first 
and  second  kind  [7].  These  are  defined  by  the  relations 


with,  in  particular,  {*}  :=  £,(-!)*©(*  -  *)*//!-  Finally,  let  Ar  be  the  forward  difference 
operator  in  r,  i.e.  A rfT  —  fr+ \  —  fr. 

Remark.  Note  our  use  of  t  versus  l.  When  referring  to  a  fixed  and  specified  agent,  i  is 
used  as  a  designator,  but  l  is  the  generic  subscript.  Hence  £  should  be  interpreted  as  “the 
specific  fixed  agent  a/’. 


The  following  acronyms  are  used  throughout: 


LTP  -  Law  of  Total  Probability  ([4]); 

CL  -  Combinatorial  Lemma  (section  B.2); 

PL  -  Probability  Lemma  (section  2.1); 

EDCI  -  Exponential  Decay  of  Common  Interest  (section  1.1) 


1.2  Assumptions  and  problem  statement 

The  derivation  of  the  crossing  network  exchange  model  is  based  upon  the  unbiased 
assumption  which  holds  throughout  this  paper  and  is  comprised  of  the  following  two 
parts: 
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(Al)  For  every  agent  o/,  each  resource  name  in  the  pool  Pi  is  equally  likely  to  appear 
in  their  state  list; 

(A2)  For  every  resource  name,  Pr(s,-  =  +1)  =  Pr(si  =  —1)  = 

The  unbiased  exchange  problem  is  to  find  the  distribution  of  J  (respectively,  Je,  Jis)  under 
the  assumptions  (Al)  and  (A2).  This  question  is  not  simply  answered  for  arbitrary  L  so 
the  focus  here  is  on  the  simpler  question  of  finding  the  means,  and  where  possible,  the 
higher  moments.  Hence,  the  following  questions  are  posed  and  answered: 

(a)  find  fi  :=  E[J}- 

(b)  find  pig  :=  E[Je]\  and 

(c)  find  Hh  ■=  E[Jls\ 

These  expected  exchange  rates  measure  respectively,  the  ‘system  performance’  which  is 
of  interest  for  network  management,  the  individual  agent’s  performance  and  the  perfor¬ 
mance  of  an  arbitrary  subset  of  the  agents  (a  syndicate).  They  establish  expected  unbiased 
benchmark  performances  for  systems  in  which  like-for-like  exchange  takes  place.  In  par¬ 
ticular,  for  applications  where  exchange  is  promoted  or  controlled  by  system  processes, 
this  benchmark  represents  the  ‘random  outcome’  performance  level  which  should  always 
be  exceeded.  If  such  performance  is  not  reached,  both  agents  and  system  management  or 
regulators  would  have  cause  for  concern. 
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2  Preliminary  results 

This  section  gives  some  elementary  discrete  probability  results  which  follow  from  the 
unbiased  assumptions  above.  The  first  result,  called  the  Probability  Lemma  (PL),  is 
pivotal  for  all  that  follows  in  that  it  facilitates  the  use  of  the  results  of  Deming  and 
Glasser  [1],  These  results  are  surveyed  in  section  2.2. 


2.1  A  Probability  Lemma 

The  equi-probable  assumption  (A2)  implies  a  useful  anonymity  when  dealing  with  state 
lists.  The  results  here  use  this  assumption  to  relate  the  number  of  exchangeable  names 
between  state  lists  to  the  number  of  common  names  among  those  state  lists.  Given  (A2), 
the  link  is  a  simple  counting  exercise. 

Let  Ut  (respectively,  Us)  be  the  random  variable  of  the  number  of  names  common  to 
the  t  (respectively,  s)  given  state  lists  of  agents  a*,,---  ,ait  (respectively,  ahl , •  •  •  ,ahs) 
and  let  Jt  (respectively,  Js)  be  the  random  variable  of  the  number  of  exchangeable  names 
among  these  Ut  (respectively,  Us)  common  ones.  The  abbreviations  P{jt\ut)  =  Pr (Jt  = 
jt\Ut  =  ut),  P(ut)  =  Pr (Ut  -  ut),  etc  are  used  to  simplify  the  notation. 

Probability  Lemma  (PL).  Under  (A2),  the  following  hold: 

(i)  If  a  given  t  state  lists  have  ut  names  in  common,  then 

P(jt\ut)  =  ( 2*-1  -  1)A2-P-1K 

(ii)  The  generating  functions  Jt(z)  :=  P{jt)zjt  and  Ut{z)  :=  J2U,  p(ut)  zUi  are  re¬ 
lated  by 


Jt(z)  =  Ut(2i-t  +  (l-21~t)z ), 

so  that 

El[Jt]k]  =  (l-2'-t)kE[[Ut}k],  k  =  0,1,2,..., 
E[Jt}  =  (l-2l~t)E[Ut], 

E[J? }  =  (1  -  2l-t?E[U?}  +  21_t(l  -  2l~t)E[Ut}- 


(in)  For  i\  ■  ■  ■  it  ^  h\  ■  ■  ■  hs  as  combinations,  the  covariance  of  Jt  and  Js  satisfies 

E[JtJs\  =  (1  -  21~~t)(l  -  2 l~s)E[UtUs\. 

Proof.  The  proofs  of  (ii)  and  (iii)  rely  on  (i)  which  is  proved  first.  Suppose  that  t  lists  have 
ut  names  in  common,  then  under  the  equi-probable  assumption  (A2),  P(jt \ut)  is  just  the 
number  of  signed  configurations  of  the  t  lists  with  ut  common  names  for  which  jt  names 
are  exchangeable  divided  by  the  total  number  of  possible  signed  configurations  of  the  lists. 
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We  need  only  consider  the  configurations  among  the  common  names  in  computing  this 
ratio.  Let  $1, . . .  ,$,i(  be  the  common  names  and  consider  the  configurations  of  the  lists: 

{s?# i,... 

There  are  2tUt  such  possible  configurations  for  6  {+1,-1}.  Suppose  there  are  jt 
exchangeable  names.  There  are  exactly  ways  of  choosing  the  exchangeable  names  each 
of  which  has  23t  possible  ways  of  choosing  the  Si’s  for  these  jt  names  and  2Ut~3t  possible 
ways  of  choosing  the  .s2’s  for  the  remaining  ut  —  jt  non-exchangeable  names.  Hence, 
2“‘  3t  is  the  number  of  possible  configurations  of  one  list  with  jt  exchangeable 

names.  Now  fix  a  configuration  of  the  first  list  and  consider  the  remaining  t  —  1  lists. 
For  each  of  the  jt  names,  ,dl  say,  that  are  exchangeable,  there  are  2l  ~ 1  —  1  configurations 
of  the  other  t  —  1  lists  in  which  that  name  has  at  least  one  list  with  an  opposing  sign 
si  *n  some  list,  i.e.  2 t  1  possible  configurations  less  the  one  configuration  where  has 
the  same  sign  in  all  the  remaining  t  —  1  lists.  Hence  there  are  (2*_1  —  l)-7*  possible 
configurations  in  which  jt  names  are  exchangeable  with  the  first  list  fixed.  The  number 
of  possible  configurations  with  jt  exchangeable  names  among  the  ut  common  names  is 
therefore  ujt  :=  (2(  1  —  1  )Jt2“f(“‘),  and  since  —  2tut ,  then  P(jt\ut)  is  the  ratio  of 

these  two  numbers  as  required  by  (i). 

To  derive  (ii),  apply  the  LTP  to  write  the  generating  function  for  Jt  as 

Mz)  =  z3tp(. itWt)P{ut). 

jt,Ut 

Now  use  the  explicit  form  in  (i)  to  write  z^P{jt\ut)  =  {2l~l  +  (1  -  21~t)z)Ut  so  that 
Jt{z)  =  P(ut){21-t  +  (1  -  21~t)z)Ui  =  Ut{ 21_t  +  (1  -  21~t)z). 

Ut 


Since  Jt{l  +  z)  generates  the  factorial  moments  of  Jt,  and  Jt{  1  +  z)  =  Ut{\  +  (1  -  2 
then  the  factorial  moments  are  related  by  scaling  with  powers  of  1  —  21~t.  The  cases  k  =  1 
and  k  —  2  then  imply  the  given  relations  between  the  first  and  second  moments. 

To  prove  (iii),  again  apply  the  LTP  and  Bayes’  rule  to  get 


P(jt\js)  ^YP^\Ut)P(Ut\js)  =  P(jt\ut)P(ut\us)P(us\js) 

Ut  Ut,Us 


Y  p(p\ut ) 


P{ut,us)  P(js\us)P(us) 

P(Us)  P(js ) 


and  so 


p(jt,js)  =  Y  P(it\ut)P{js\us)P{ut,us). 

Ut,Us 
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Now  taking  the  required  expectation  gives 

E[JtJs\  =  EE  3tP{jt\ut)jsP{js\us)P{ut,us) 

Ut,us  jt,js 

=  (1  -  21  ()(1  -  21~s)  utuaP(ut,us), 

Ut,Us 

where  again,  the  expression  in  (i)  has  been  used  to  evaluate  the  j-sums.  □ 

Observe  that  the  conditional  probability  in  (i)  of  the  PL  is  in  fact  a  Binomial  distribution 
on  ut  trials  with  the  trial  success  probability  p  =  1  -  21-*  and  the  trial  failure  probability 

1  —  p  =  21-<. 


2.2  The  Deming-Glasser  distribution 


Before  the  advent  of  modern  computing  resources  for  record  keeping  and  manipula¬ 
tion,  sampling  methods  to  assist  with  counting  problems  were  a  topical  issue  of  research. 
For  example,  the  works  of  Deming  and  Glasser  [1]  and  Goodman  [2]  were  motivated  by 
the  problem  of  comparing  large  lists  of  names  to  determine  the  size  of  common  member¬ 
ship  through  the  comparison  of  small  random  samples  taken  from  those  large  lists.  Two 
particular  motivating  applications  for  this  work  were: 


•  determining  the  number  of  names  common  to  different  lists  of  social  welfare  recipi¬ 
ents; 


•  determining  the  number  of  books  common  to  different  library  catalogues. 


In  their  day,  these  posed  significant  problems  because  of  the  lack  of  ease  of  comparing 
large  collections  of  paper  records.  This  section  gives  a  brief  summary  of  those  of  Deming 
and  Glasser’s  results  from  [1]  which  are  used  below  in  the  crossing  network  model. 

Consider  t  overlapping  pools  of  names  Pi, . . .  ,Pt  of  respective  sizes  Nu  . . .  ,Nt  with 
C  names  in  common  (see  Figure  1).  For  t  independent  respective  random  samples  of 
sizes  n i,  ■ . .  ,nt  taken  from  these  pools,  Deming  and  Glasser  [1]  showed  that  the  random 
variable  c  of  the  number  of  names  common  to  the  t  samples  has  the  distribution 


fi) 

ii  fit,\ 

1=1  \k  ) 


C  =  0,1,. 


,  mm(n;,  C ) 


which  is  called  the  DG  distribution.  Some  simple  rearranging  of  terms  shows  that  the 
generating  function  Pt(z)  :=  Y^cPt(C\c)zc  is 


Pt{z)  =t+iFt{-C;-nu...  ,  -nf;  -Nu  . . .  ,  —Nt;  1  —  z), 

where  pFq  is  the  generalised  hypergeometric  function  [5].  This  identifies  the  DG-distribution 
as  a  generalised  hypergeometric  factorial  moment  distribution  (ghfd)  [3].  Since  Pt(l  +  z) 
generates  the  factorial  moments  :=  £?[[c]*],  then  immediately, 

R[fc]  =  iC\k  II  k  =  0, 1, ,  min (m,  C ). 
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Alternatively,  since  the  factorial  moments  can  be  identified  by  applying 

the  Lemma  of  section  B.2  with  bk  =  pj k]/h\  and  pr  =  Pt{C ;  r).  The  derivation  of  Pt{C;  c)  is 
a  straight  forward  application  of  the  LTP  (cf.  section  3.1).  Consequently,  Pt(C\  c)  satisfies 
a  recursion  which  is  now  discussed. 


If  an  additional  ( t  +  l)-st  independent  random  sample  of  size  rp+i  is  taken  from  a  pool 
-Pf+i  °f  size  Nt+ 1,  still  with  C  names  in  common  with  the  pools  Pi, . . .  ,Pj,  then,  by  the 
LTP, 


C  fk\fNt+1-k\ 

pt+\{C]c)  -  y:  —  C  Pt{C;k),  c  =  0,1,...  , min(n/, C). 

k=0  Vnt+ii 


(1) 


Moreover,  for  each  r  -  0, . . .  ,  t ,  the  recursion 


c 

Pt+I  (C;c)  =  Y^  pr+ 1  (k\  c)Pt-r  (C;  k ) 

k—0 

follows  from  the  implied  repeated  sum  representation  of  Pi+^Cjc)  derived  by  iterating 
(1).  In  this  recursion,  P0{C-,k )  =  8kC  and  Pi(C;fc)  =  (J)  More  explicitly, 

augmenting  the  notation  to  include  all  of  the  omitted  parameters, 


0  Q  ^  Pi  shows  this  last  recursion  in  detail.  The  main  property  of  the  DG-distribution 
required  here  is  its  mean: 

*H-4]=cn$.  <2) 


2.3  A  special  exchange  problem 

This  section  considers  a  more  restrictive  version  of  the  resource  exchange  problem  in 
order  to  illustrate  the  use  of  the  DG-distribution  in  obtaining  explicit  results  for  the  number 
of  exchangeable  names  among  agents.  The  special  problem  solved  here  is  to  determine 
the  distribution  of  the  number  of  exchangeable  names  only  within  those  names  common 
to  all  of  a  given  set  of  resource  pools. 

Let  Jt  be  the  random  variable  of  the  number  of  exchangeable  names  in  the  state  lists  of 
agents  ai, . . .  ,  at,  2  <  t  <  L,  among  those  C  (drop  the  indexing  notation  C\...t  of  section 
1.1)  names  common  to  the  pools  P\ , . . .  ,  Pt  and  let  c  be  the  random  variable  of  the  number 
of  names  common  to  all  t  state  lists.  Then  by  the  LTP  and  the  PL, 

Pr  (Jt=j)  =  Y/P(j\c)Pt(C;c) 

C 

=  (2t_1-i  )JX(c)2“(t-1)c^(tf;c), 

c= 0  ''IJ 
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Figure  1.  A  schematic  for  t  —  3  showing  three  samples  of  sizes  ni  taken  from  overlapping 
pools  of  size  Nt  -  dark  shading  is  the  random  variable  c  while  light  shading  indicates  C. 


where  Pt{C\c )  is  the  DG-distribution  with  the  parameters  for  agents  ai, . . .  ,a.t.  Forming 
the  sum  JT  z^Pr(Jt  =  j),  and  evaluating  the  j-sum,  it  follows  that  the  generating  function 
for  Jt  is 

Jt{z)  =  ,+i  Ft(-C;  -nu  ...  ,  -nt;  -Nu. . .  ,  -Nt;  i(l  -  z)) 

ajid  the  central  moments  fi^  oi  Jt,  derived  from  the  factorial  moment  generating  function 
Jt{  1  +  z)  [4,  3],  are 


with,  in  particular,  the  mean  and  variance  being  given  by 


(-£)',  k  —  2,3, 


p=(i-2-)cn^ 

i=i  1 


a2  =  fi2  =  A[1  -  id  +  (1  -  21~~t)(C  -  1)  JJ 


i=i 


nt  -  1  - 

N,  -  1- 


Thus,  the  distribution  and  the  moments  of  the  number  of  exchangeable  names  in  those 
names  common  to  all  pools  of  a  given  set  of  agents  can  be  written  down  explicitly  using 
the  DG-distribution.  The  difficulty  in  determining  the  total  number  of  exchangeable  names 
overall  lies  in  not  double  counting  these  names  when  they  lie  in  multiple  overlaps  of  the 
pools  of  many  agents  (cf.  Figure  1). 
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Finally,  it  follows  from  the  expressions  for  p  and  cr2  that  for  an  observation  j  of  Jj,  the 
estimator  C  j{  1  —  21  f)  1  A;  1  is  an  unbiased  estimator  of  C  with  standard  error 


var  C  = 


1— 1  r 


Ih 


l  +  (C'-l)(l-21-i)J] 


ni 


Ni-  1J 


-C2. 


and  an  unbiased  estimator  of  this  standard  error  is: 


V  :=  C 


(1 


n 


N,  -  1 

ni  -  1 


+  C2 


II A' 


Nt-1- 
ni  -  1  - 


(ef.  [4]).  Clearly,  the  agents  Oi,  ■  •  •  ,  at  can  be  replaced  by  any  subset  ap , . . .  ,  alt  of  the  L 
agents  and  the  above  results  still  hold  with  the  replacements  A,  Aiz  and  C  ->  Ch...it. 
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3  A  few  agents 

We  begin  a  study  of  the  unbiased  exchange  problem  by  considering  the  two  cases  L  =  2 
and  L  =  3.  Apart  from  introducing  and  illustrating  the  techniques  used  in  the  general 
case,  these  cases  are  more  easily  solved  to  reveal  more  moments  of  the  exchange  random 
variables  J  and  Jt.  In  particular,  the  basic  case  L  =  2  is  completely  solvable  through  the 
characterisation  of  the  exchange  distribution  for  J,  while  for  the  L  =  3  case,  only  certain 
common  name  joint  distributions  and  the  first  and  second  moments  are  calculated. 


3.1  The  2  agent  case 

The  two  agent  case  is  fundamental  to  exchange  in  many  agent  scenarios:  any  exchange 
ultimately  takes  place  between  pairs  of  agents  as  they  are  compared  or  meet  in  a  network. 
This  case,  therefore,  forms  an  intuitive  building  block  of  more  complicated  situations.  For 
L  =  2,  the  questions  (a),  (b)  and  (c)  of  section  1.2  collapse  to  the  one  question  since  every 
exchangeable  name  in  one  list  is  in  the  other  (i.e.  J  =  Jt).  Since  there  is  only  one  overlap 
of  resource  pools,  the  solution  of  the  exchange  problem  for  this  case  has  already  been  given 
in  section  2.3  (set  t  =  2).  However,  a  separate  derivation  is  instructive.  After  computing 
the  exchange  distribution  and  all  of  its  moments,  the  questions  of  small  mean  behaviour 
and  the  estimation  of  the  number  of  names  common  to  the  two  resource  pools  are  briefly 
considered. 


3.1.1  The  exchange  distribution  and  its  moments 


To  aid  in  simplification,  the  notation  here  departs  slightly  from  that  in  section  1.1. 
Consider  two  agents  a\  and  02  with  state  lists  of  length  n  and  m  respectively  consisting  of 
respective  samples  of  resource  names  from  the  pools  Px  and  P2  which  are  of  respective  sizes 
N  and  M.  The  quantities  7r,-(l,  2)  are  written  simply  as  7 q  and  Xx  =  n/N  and  A2  =  m/M. 

Let  j  be  the  number  of  exchangeable  names  and  suppose  that  Px  and  P2  have  C  names 
in  common.  Write  pj  =  Pr(J  =  j).  In  this  special  case,  a  complete  characterisation  of  pj 
is  readily  available.  In  fact,  we  show  that  for  two  agents, 


Pj  = 


j  =  0, . . .  ,  min(m,  n,  C). 


The  proof  of  the  formula  for  pj  is  an  application  of  the  LTP.  Define  the  events: 


Bc  =  n-list  contains  c  of  the  C  common  names,  0  <  c  <  C] 

Ak  =  m-list  has  k  common  names  with  the  n-list,  0  <k  <m\ 

Dj  =  j  exchanges  occur  between  the  n-list  and  m-list,  0  <  j  <  m. 
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Then  using 


c 

Pv(Ak)  =  J2Pr(Ak\Bc)Pv(Bc), 

c— 0 
m 

Pr(^)  =  EPr(^i^)pr(^). 

k= 0 

with  Pr(Bc)  =  ©_(";f)/("),  Pr(4t|Bc)  =  ©("_!)/(")  and  Pv(D,\Ak)  =  2“‘(‘)  (* 

successful  probability  |  Bernoulli  trials,  cf.  t  =  2  in  the  PL),  gives  the  expression  for 
Pj.  The  hypergeometric  probabilities  appear  as  a  reflection  of  the  (equally  likely)  selection 
without  replacement  of  resource  names  from  the  pool.  The  distribution  P2(C;  k)  :=  Pr(i4*), 
which  gives  the  number  of  common  names  between  the  two  lists,  is  just  the  DG-distribution 
for  two  lists  (see  section  2.2). 

The  distribution  {pj}  is  a  generalised  hypergeometric  factorial  moment  distribution 
(ghfd)  [3]  and  is  a  (probability  |)  Bernoulli  generalisation  of  the  common  name  distri¬ 
bution  P2(C]k)  in  [1].  Note  the  symmetry  under  simultaneous  interchange  of  m  with  n 
and  M  with  N.  The  distribution’s  properties  are  summarised  through  considering  the 
generating  function  P{z)  :=  YljPjzP  Let  pk  :=  E[(J  —  p)k]  be  the  central  moments  with 
p  =  P[J]  the  mean  and  a 2  the  variance.  The  generating  function  is 

P{z)  =  3F2(-C,  —n,  -m;  —N,  —M;  |(1  -  z)), 

and  the  central  moments  pk  are  given  by 


Pk 


E2_i 


Wi[n]i[C]j 

[M]i[N]i 


E 

1=0 


k  =  2,3,...  . 


In  paiticular,  the  mean  (p),  variance  (<72)  and  (Fisher)  skewness  (7  :=  ps/a^)  are  given 
by  . . 


p  —  |AiA 2C, 

a2  =  p(l  - /i  +  |(C  -  1)tti),  (3) 

1  [1  —  3//  +  2p2  +  |(1  —  p)[C  —  1)717  +  \{C  —  1)(C  —  2)7t2] 

\[p  [1  —  /Z  A(C  -  l)7T1]3/2  • 

Although  seemingly  complicated,  these  expressions  for  the  moments  follow  easily  from  the 
generating  function.  Since  P(1  +  z)  generates  the  factorial  moments,  then  identifying  the 
coefficients  of  zk  in  P(l  +  z)  shows  that  p[k]  ■=  E[[J]k ]  =  2-k[m}k[n}k[C}k/{[M]k[N]k)  from 
which  standard  moment  relationships  [3,  4]  give  the  expressions  for  the  central  moments 
Pk- 

Observe  that  all  the  results  here  display  the  symmetry  identified  earlier  expressing  the 
freedom  to  re-label  the  agents.  The  distribution  {pj}  also  has  a  number  of  reductions 
which  can  be  identified  from  the  corresponding  P(z)  -  these  are: 


•  p  C  —  M,  then  P(z)  —  2Pi(— n,  — m;  —N;  i(l  —  z))  which  characterises  the  hyperge¬ 
ometric  distribution  (with  parameters  n,  m  and  N)  generalised  by  the  (probability 
5)  Bernoulli  distribution  (similarly  for  C  =  N); 
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•  If  n  —  TV,  then  P(z)  —  2-^1  (~C,  —  m;  —  M;  |(1  —  z))  which  is  the  hypergeometric  dis¬ 
tribution  (with  parameters  C,  m  and  M)  generalised  by  the  (probability  1)  Bernoulli 
distribution  (similarly  for  m  =  M); 

•  If  n  —  N  and  m  =  M,  then  P(z)  —  \Fq(  — C;  |(1  —  z))  which  is  a  binomial  distribution 
on  C  trials  with  success  probability 

•  If  m  —  M  and  C  —  TV,  then  P(z )  =  \Fo{— n;  ±(1  —  z))  which  is  a  binomial  distribution 
on  n  trials  with  success  probability  i  (similarly  for  n  =  TV  and  C  =  M). 


j  n 

(a)  (b) 

Figure  2:  For  the  parameter  choices  m  —  n,  N  —  M  =  100  and  C  =  50,  (a)  plots  the 
distributions  {pj}  for  the  three  values  n  =  30,  50,  100  (left  to  right)  while  (b)  plots  the 
moments  p,  a 2  and  7  as  functions  of  n. 

Figure  2  shows  example  distributions  {pj}  in  (a)  and  the  first  three  moments  in  (b). 
By  fixing  all  but  one  of  the  size  parameters  (nj,  one-dimensional  plots  result.  In  these 
examples,  the  two  agents  are  similar  in  that  m  —  n  and  M  —  TV.  Figure  2(a)  illustrates 
the  transition  to  the  binomial  distribution  as  n  increases,  while  (b)  shows  the  quadratic 
(in  n)  mean,  an  inflection  in  the  variance  and  the  rapidly  decaying  skewness. 


3.1.2  Small  mean  behaviour 


The  mean  and  variance  approximately  agree  for  small  p.  This  is  seen  by  writing 


9  9 

a  =  p  +  pl 


(1  —  m  1)(1 


•  n 


-1 


L(1  -M-!)(l  -  TV-1) 


(i  -c- 


which  also  illustrates  that,  for  p  small  (i.e.  if  any  of  n/TV,  m/M  or  C  are  small),  the  level 
curves  mn  =  const  of  /x  also  approximately  preserve  the  variance.  The  degree  of  agreement 
of  the  two  curves  is  shown  in  Figure  2(b). 

The  exchange  distribution  exhibits  a  favourable  positive  skewness  which  increases  for 
decreasing  mean  p  (see  Figure  2(b)).  For  the  skewness,  . 


7  = 


1 


[1  +  0(/i)], 


as  p  -7  0. 
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Hence,  the  exchange  distribution  { pj }  is  highly  positively  skewed  only  for  small  mean  and 
variance.  As  fj,  increases,  7  decreases  with  in  fact  7  =  0,  when,  according  to  the  reductions 
above,  either  both  state  lists  are  of  maximum  size  (n/N  =  m/M  =  1)  or  there  is  maximal 
resource  pool  commonality  and  one  list  is  of  maximum  size  (e.g.  n/N  =  1  if  C  =  M).  This 
shows  that  as  state  list  size  increases,  favourably  skewed  exchange  outcomes  gradually 
disappear.  The  rapid  rate  of  this  decrease  is  illustrated  in  Figure  2(b). 

3.1.3  Estimation  of  common  names 

Estimating  the  number  of  common  names  (C)  to  the  resource  pools  from  exchange 
outcomes  is  of  interest.  If  jk,  k  =  l,--  ,  K,  are  K  independent  exchange  outcomes  for 
fixed  parameters,  then 


K 


C  = 


Kir0 


k=l 


(recall  -kq  =  A1A2)  is  an  unbiased  estimator  of  C  and,  using  (3),  has  a  standard  error 


var  C  = 


7T0  K 


1  + 


{C  -  l)7Tl 


K ' 


Moreover,  an  unbiased  estimator  of  this  standard  error  is 


V 


c_ _ (2  -  it  1 ) 

7T0  (K  +  TTi/tTo  -  1) 


C_ 

K 


(2/tto  -  1) 


if  7ri/7r0 


_ a _ \ 

K  —  1  +  7T !  /  7T0  / 


These  results  area  a  AT-sample  version  of  the  t  =  2  estimators  in  section  2.3.  Estimates  of 
other  parameters  can  be  made  similarly. 


3.2  The  3  agent  case 


Consider  three  agents  a\,  02,  03  whose  resource  pools  Pi,  P2,  P3,  are  described  by 
overlap  parameters  C/m  for  pairs  /,  m ,  and  have  C  names  in  common  to  all  three  pools. 
The  indices  l ,  m  and  q  are  used  to  denote  different  agents.  The  L  =  3  case,  being  more 
complex  than  L  —  2,  serves  as  an  introduction  to  the  method  used  for  the  general  L 
case  in  section  4,  however,  more  is  explicitly  calculable  than  for  the  general  case.  In 
particular,  some  joint  common  name  distributions  can  be  derived  which  in  turn  allow  the 
easy  computation  of  the  variances  of  the  exchange  random  variables  J  and 

3.2.1  Common  name  distributions 
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We  begin  by  deriving  the  joint  distribution  of  the  number  of  common  names  between 
two  and  three  state  lists.  Let  Qm  be  the  random  variable  of  the  number  of  names  common 
to  the  state  lists  of  a;  and  am,  c  the  number  of  names  common  to  all  three  state  lists  and 
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dim  ■=  Cim  -  C  the  number  of  names  uniquely  common  to  the  state  lists  of  a,  and  am.  Let 
r  be  the  number  of  names  among  the  C  names  common  to  all  pools  which  are  common  to 
the  lists  of  a;  and  am,  then  random  choice  without  replacement  implies  that 


P{r\cim)  = 


(c'rm)  (Cl  m  Qm  ^ 


and  so  P(r,  c/m) 


Im  C[  rn  ^ 


P2(Q  m  i  C/m)* 


Consequently,  application  of  the  LTP  in  the  form  P(c|c/Tn)  =  £r  P(c|r)P(r|c,m)  gives  that 


P(c\cim)  =  Y 


0  (nq-c )  (%m)  ( 


Clm  Clr 
C—r 


r= 0 

m 


(  q) 

\TlqJ 

(Na- 


C Ccm ) 


P{pi^lrn)  —  ^  ] 


(c)  (rig—c)  (Crm)  (C!c-rm) 


r= 0 


r) 

\na  ) 


{%-) 


Qm)- 


To  check  that  P(c,  clm)  has  the  correct  reductions,  £c  P(c,  c/m)  =  P2(Om;  clm )  is  obvious, 
however,  summing  over  c/m  requires  the  identity  (for  C  <  C/m) 

Y  ■  — /C-JV  f-AQm-.Cm)  =  P2(C;r) 

C(m  VC/ 


to  hold.  This  is  just  an  application  of  the  dg  recursion  (1)  in  which  the  parameter  Ctm 
cancels  and  reduces  the  order  from  t  =  3  to  t  =  2.  The  remaining  sum  over  r  then  produces 
P3 (C]c)  as  per  (1)  so  that  X^cim  P{ci  cim)  =  Ps(C;c).  Prom  P(c,  C/m),  we  compute  that 
E[cclm]  =  \q(C/Clm)E[cfm},  and  hence, 


E[cClm]  —  A;AmA9C[  1  +  (C;m  —  l)7Ti(/,m)]  (4) 

follows  since  P[cfm]  =  A/AmC/m[l  +  (C/m  —  l)7Ti(Z,m)]  is  the  second  moment  of  the  DG- 
distribution  P2(C/m;Qm)  (see  section  2.2). 

The  joint  distribution  of  C[m  and  c;g  can  also  be  found  by  the  application  of  the  LTP 
and  Bayes’  rule.  These  details  are  relegated  to  an  Appendix  (section  B.l),  however,  the 
result  is: 


P(c/m)  clq) 


(k\(Nq-k 

E\Clq/  \Tlq- 

('M 

\na  / 


q  clq'  \C[ 


)(„')( 


Nm~ 
nm  ~clm 


71/  —  S 


fc,r,s 


(/Vm) 

\nm) 


o 


o 


(r)  feD  Ofe*)  (n,0 

fi-)  (t)  B(M‘ 


(5) 


This  distribution  too  satisfies  the  appropriate  reductions  such  as  £  P(ctm,  cl9)  =  P2(C';9;  q9), 
and,  when  C  =  0,  it  factors  as  P2{Cim ;  c;m)P2(Qg;  ciq)  to  reflect  the  resulting  independence 
of  c(m  and  q,.  The  only  need  for  P(c/m,cf,)  is  in  the  computation  of  P[c,mc,9]  which  is 
required  below.  Evaluation  of  this  expectation  requires  the  two  sums  Tlm(r)  and  Tiq(r) 
where 


Tim(r)  :=YkP 

k 


P=1,B  =  Cim. 
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The  Appendix  (section  A.l)  gives  a  method  for  evaluating  this  sum  for  p  an  arbitrary 
non-negative  integer  and  arbitrary  B  >  C.  The  result  for  p  =  1  and  B  =  Cim  (see  section 
A. 2)  is: 


Tlm(r) 


(Qm  z  ££ V  +  W  -  clm)r  0(";_3 


N,-C  (";)  ' 

This  result  is  the  key  part  of  the  computation  of  E\cimciq\  which  then  turns  out  to  be: 

P\clmclq]  —  ^Am\|c[l  +  (C  —  V)tX\  (Z)]  +  C7Ti(/)(C';m  -f-  Clq  —  2 C) 

1  -  TT\  (/) 


+  (C|m  -  C)(Clq  -  C)  7Ti(Z)  + 


(6) 


Nt-C 

Since  the  joint  distribution  P(c,  Qm)  is  known,  the  distribution  for  the  number  of  uniquely 
common  names  to  two  state  lists;  namely,  dim  :=  c/m  —  c,  can  be  written  as 


C  d  f r+c\(N-r-c\  (c+d\  (Cim—c-d\ 

V  c  /  V  7i — c  /  \r+c/  V  C—r—c  ) 


P  x{dlm  —  d)  —  P2(Cim]  C  +  d) 


c=0 


r=0 


o 


(T) 


with  the  corresponding  generating  function 


C  Clm—C+c  fr\  / N-r\ 

Dim(z)  :=^Tz-c  ^  c  ,n~c  Pim{r;z),  where 

c=0  r=c 

Clm-C+C 


Plm(r;  z)  =  J2  P2(Clm-,d) 


o 

0(Clcm--rd)  _d 


d=r 


C C  ) 


z . 


3.2.2  Combinatorial  identities 


Before  leaving  the  common  name  joint  distributions,  it  is  interesting  to  note  some 
combinatorial  identities  which  follow  from  the  LTP  written  in  the  forms 

P(c\ °lq)  =  yy  P(c\cim)P(cim\ciq)  (7) 

^Im 


and 


P(cim\Clq)  =  'SjTjP{cim\c)P{c\ciq). 

C 

When  Cim  =  C/9  =  C,  P(cim,cig )  and  P(c,  c/m)  respectively  reduce  to 


c 

P{clmi  clq)  —  ^  ^ 
k= o 


(  k  w  Nm-k  \  (C\  (N,-C\  (M/M) 

\7lffi  C/m/  \k )  \7li — k )  'C/q/  \TLq  Clq  / 


and 


(8) 


Qm  )  — 


(c,r)  (X~cl) 


p2  (C;  Cim). 
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For  this  special  case,  application  of  the  forms  (7)  and  (8)  produce,  respectively, 


where  P2(C;r)  :=  P2(C;n,,nm;iVj,./Vm;r)  and  P2(C;s)  :=  P2{C]nl,nq- Nh  Nq;  s).  Revert¬ 
ing  to  the  general  case  C;m  ^  Ciq  7^  C,  then  applying  (8)  and  using  the  last  identity  above 
to  replace  the  summation  over  c  produces 


P{clmiclq)  PiiGlm'i  Ciy 

(Crl)  (Q  m 

r  (4) 


E 

r,s,k 


i)Pi  (Clq>  Ctq) 

(*)(£:*)  ©(©T)  (IKK)  ©ff;T 

:)  ©•)  Pi(c;s)C:) 


P2(C',r)( 


Nn 

Pm  > 


© 


It  remains  to  investigate  the  novelty  of  these  identities  in  terms  of  their  place  within 
existing  frameworks  [6,  7]  . 


3.2.3  The  exchange  problem 

If  j  (respectively,  ji)  is  the  number  of  exchangeable  names  among  the  state  lists  of  all 
three  agents  (respectively,  in  the  state  list  of  agent  af),  write 

j  ~  ^2  jlm  +  Jl23) 
l,m 

3i  ~  Jim  +  jig  +  j  123, 

where  jim  is  the  number  of  exchangeable  names  among  the  d/m  names  uniquely  common 
to  the  state  lists  of  agents  a/  and  am  and  j\23  is  the  number  of  exchangeable  names  among 
those  names  common  to  all  three  state  lists.  The  decomposition  into  exchange  among 
uniquely  common  names  conveniently  separates  common  names  into  disjoint  classes  of 
maximum  overlap  among  lists.  The  same  principle  is  used  later  in  the  general  case.  From 
(ii)  of  the  PL, 


E[j\  =  (1  -  2  ')  Y,  E[dlm]  +  (1  -  2~2)P[C], 

l,m 

E[je\  -  (1  -  2-l)E[dtm]  +  (1  -  2~1)E[deq\  +  (1  -  2“2)£;[c]. 

Since  dim  =  C;m  —c  and  E[cim}  =  A/AmC(m,  E[c]  =  A1A2A3C,  then  the  means  n  =  E\j]  and 
\H  =  E[j(]  follow  as 


R  I  ^  .  A/ AmCjm  |Ai  A2A3C 

l,m 

He  =  |AfAmC£m  +  ^i^qC(q  —  IA1A2A3C. 
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Comparing  these  means  with  (3),  the  expected  outcome  is  a  sum  of  two-agent  means  less  1 

a  correction  for  mutual  overlap. 

To  find  the  corresponding  variances  a 2  and  cr2,  E[j2]  and  E[j2]  are  computed  using 
the  common  name  joint  distributions  from  section  3.2.1.  Expanding  the  corresponding 
expressions  gives 


3  XX  ^  'y  '  jlmjlq  T  2ji23  y  jlm  T  ji23i 
l,™  l^m,q  l,m 

3e  ~  3im  +  dlq  +  2 jlmjlq  +  % j  123 jtm  +  % j  123 jiq  +  ji23, 

from  which,  by  taking  expectations,  it  follows  from  (ii)  of  the  PL  that 

E[j2]  =  (1  -  2-1)2  ]T  E[dfm]  +  2“1(1  -  2-1)  E[dlm)  +  2(1  -  2~1)(1  -  2~x)  £  E[dlmdlq] 

(,m  l^rn,q 

+  2(1  -  2-2)(l  -  2-1)  Y,  E[cdlm]  +  (1  -  2-2)2£[c2]  +  2“2(1  -  2~2)E[c] 

l,m 


and 

E\j]}  =  (1  -  2-1)2£[4n]  +  2~1(1  -  2-1)E[dem]  +  (1  -  2-1)2E[djg}  +  2~1(1  -  2 ~l)E[dtq\ 

+  2(1  -  2— 1 ) (1  -  2 -l)E[dimdiq}  +  2(1  -  2_2)(1  -  2-1)E[cdim] 

+  2(1  -  2"2)(1  -  2~1)E[cdtq\  +  (1  -  2“2)2J5[c2]  +  2~2(1  -  2-2)E[c\. 

As  with  the  calculation  of  the  mean,  writing  dim  =  cim  —  c  expresses  E[j2]  and  E[j2]  in 
terms  of  only  the  expectations  involving  c/m  and  c  computed  in  section  3.2.1.  Substitution 
of  the  expectations  (4)  and  (6)  then  gives 

E[j  ]  —  ^  |  A; AmC)m[l  +  \{Cim  —  1  )7Ti  (l,  m)]  +  |  A1A2A 3C(C  —  1)  y^  7Tj  (/) 

l,m  1 

~  l)7r1(/,m)  +  ^A1A2A3C(C-l)7r1(l,2,3)  -  fA^A^ 

l,m 

+  |A1A2A3{c  J2  *i(l)(Cim  +  Clq-2C)+  yr  (Clm  -  C)(Clq  -  C)[7n(0  +  {l.~ 

l^rn,q  l^m,q  '  1  ' 


and 


E{fi]  =  |A£AmQm[2  +  (Qm  -  1)t r^m)]  +  iA£A9Cy2  +  (Q9  -  IJttj (£,<?)] 

-  |AiA2A3C[(C'£rra  —  1)7Ti(£,  m)  +  (C^9  —  1)7Ti(£,  g)] 

+  ^AiA2A3C(C-  1)^(1, 2, 3)  -  iA1A2A3C+iAiA2A3C(C- 1)^(0 
+  | Ai  A2A3 ^Ctt\  (£) (C(m  +  Ctq  -  2 C)  +  (Ctm  -  C)(Ceq  -  C)[ 
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Finally,  the  variances  o 2  and  a2  follow  by  subtracting  the  squares  of  the  respective  means 
H  and  H£. 
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3.2.4  Two  special  cases 

If  C  —  0,  the  random  variables  jim  are  independent  and  the  means  and  variances 
simplify  accordingly.  Writing  =  |A,AmC/m  and  ofm  for  the  variance  of  jtm,  then 

a  ~  alm  ~  Wm[l  ~  M/m  +  2  1  (C;m  —  l)7Ti  (/,  m)] 

Z,m  l,m 


and 


c=0  —  +  oiq  —  Him[  1  -  IHrn  +  2  (Qm  -  l)ni(£,m)\ 

+  Mq[  1  -  +  2_1  (C*9  -  1)71-!  (I,  q)} 


which  expresses  the  variances  as  the  sum  of  two  agent  variances  (cf.  (3)). 

Another  interesting  case  is  when  Chn  =  C,  for  all  pairs  l  and  m,  so  that  all  agents 
share  a  common  ‘core’  of  C  resource  names  among  their  pools.  In  this  case,  the  means 
and  variances  simplify  somewhat  with  the  means  reducing  to  symmetric  polynomials  in 
the  fractions  A;  (see  also  section  4.5).  A  graphical  comparison  with  the  2-agent  case  can 
be  made  through  further  setting  nx  =  n2  =  n3  =  n  and  Nx  =  N2  =  Nj,  =  N  so  that  all 
A i  are  of  equal  value.  This  describes  three  like  agents.  Then,  using  A  for  the  value  of  each 
A i,  the  means  above  become  the  cubic  polynomials 


p  =  §CA2(1  -  A/2)  and  =  C\2(l  -  A/4). 

By  further  fixing  the  values  for  C  and  N,  these  means  and  the  variances  a 2  and  a],  as  well 
as  the  corresponding  ones  for  the  L  =  2  case  (see  (3)),  all  become  polynomials  in  A.  Both 
a2  and  of  above  are  of  degree  six  while  the  L  =  2  variance  is  of  degree  four. 

Figure  3(a)  plots  n*/ C  where  /i*  is  each  of  the  L  =  2  mean  (±A 2C)  and  //,  and  /ie 
above.  As  expected,  the  single  agent  mean  (w)  exceeds  the  2-agent  mean  but  falls  short 
of  the  mean  fi  of  all  three  agents  taken  together.  The  inflection  in  the  curve  of  p/C 
demonstrates  the  onset  of  a  characteristic  ltS-shape”  for  the  system  mean  which  becomes 
more  pronounced  with  an  increasing  number  of  agents.  Section  4.5  gives  a  more  complete 
illustration  of  this  behaviour. 

Figure  3(b)  plots  of/C  where  of  is  each  of  the  L  =  2  variance  (in  (3)),  of  and  a2. 
While  for  small  A  each  variance  is  an  increasing  function  of  both  A  and  the  number  of 
agents,  the  3-agent  variances  exhibit  a  maximum  and  then  decrease  for  A  approaching  one. 
This  implies  that  for  more  than  two  agents,  there  are  ranges  of  state  list  size  in  which  the 
exchange  rate  experiences  a  greater  variance  and  therefore  larger  deviations  from  the  mean 
are  more  probable.  From  the  agent’s  viewpoint,  it  is  desirable  to  operate  in  the  reduced 
variance  regions  where  smaller  deviations  in  the  exchange  outcome  are  predicted.  The 
extent  and  location  of  the  variance  s  maximum  can  be  altered  by  changing  the  parameters 
N  and  C. 
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(a)  (b) 


Figure  3:  For  the  values  N  =  100  and  C  =  50,  (a)  plots  the  three  means  p/C,  pifC  and 
jA2  (L  —  2)  ( top  to  bottom)  while  (b)  plots  the  variances  a2 /C  and  (J2/C  as  the  curves 
exhibiting  maxima  (a2  >  a2)  and  the  variance  from  (3)  for  L  =  2,  all  as  functions  of  the 
state  list  fraction  A  =  n/N  (see  the  text). 


( 


4  The  general  case 

This  section  derives  the  mean  of  the  number  of  exchangeable  names  for  an  arbitrary 
number  L  of  agents.  The  results  also  include  the  mean  for  an  arbitrary  sub-collection 
(syndicate)  of  agents  from  the  total  of  L  agents.  In  particular,  the  case  of  all  L  agents 
considered  together  (the  system  result)  and  the  case  of  a  single  agent  (the  individual 
result)  are  special  cases  of  interest.  A  system  administrator  or  resource  controller  would 
be  interested  in  the  expected  exchange  rate  for  all  L  agents  together,  while  each  of  the 
independent  agents  is  more  concerned  with  their  own  personal  exchange  rate  as  a  way  of 
quickly  gaining  or  disposing  of  resources.  In  applications,  crossing  network  regulators  and 
user  agents  are  interested  in  achieving  maximal  exchange  rates.  The  results  derived  here 
set  minimum  performance  benchmarks  for  these  applications. 

In  the  next  three  sub-sections,  all  agents  (section  4.1),  a  single  agent  (section  4.2)  and 
an  arbitrary  subset  of  s  of  the  L  agents  (section  4.3)  are  considered.  Although  the  results 
in  section  4.3  can  be  specialised  to  produce  the  single  and  all  agents  cases,  the  derivation 
in  section  4.1  is  important  in  establishing  both  the  general  method  and  a  repeatedly  used 
expression  for  the  number  of  uniquely  common  names  to  an  arbitrary  collection  of  state 
lists  (see  (12)). 


4.1  All  agents 

Consider  now  the  problem  of  finding  p  =  E[J ]  for  the  whole  collection  of  L  agents. 
By  considering  all  agents,  evaluating  p  answers  the  question  of  a  system  performance 
benchmark  for  the  crossing  network.  The  notation  here  follows  that  of  section  1.1. 

We  will  show  that  the  expected  number  of  exchangeable  names  p  among  L  unbiased 
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\ 


I 


agents  is 

L 

H  =  E[J]  =  -  21_t)  E  Civ..it •  •  •  K-  (9) 

t= 2  *i  -  it 

The  basic  idea  of  the  derivation  is  to  partition  the  common  names  among  state  lists 
into  disjoint  sets  of  uniquely  common  names  to  A;-cohorts  of  lists  and  then  apply  the  PL. 
This  avoids  double  counting  the  common  names  which  lie  in  multiple  state  lists.  The 
total  number  of  exchangeable  names  is  then  simply  the  sum  of  the  exchangeable  uniquely 
common  names  (cf.  section  3.2.3). 

Let  Xk  •—  Sq-.-i*.  denote  the  number  of  names  uniquely  common  to  all  combi¬ 

nations  of  k  state  lists.  We  first  show  that 

**  =  D-  <10> 

t-k  ^  '  *i ---it 

and  so,  by  (2), 

E[xk]  =  Y^~^t~k(k]  Civ"itXii  Xir 

t-k  ^  '  ii  -it 

Since  the  number  of  names  uniquely  common  to  the  state  lists  of  the  agents  aq , . . .  ,  alk 
is  the  number  of  common  names  minus  the  number  of  uniquely  common  names  to  these 
and  all  other  state  lists,  then 

L-k 

dil—ik  =  ch  —h  ~  dix  -ikj\  -jf 

t= i  jvjip 
h>—  M 

Summing  over  all  combinations  i\  ■  ■  ■  ik  from  the  L  agents  then  gives  the  recursion 

xk=Y  ch-h  -  E  C '+kk)xt+k ’ 

h-ik  t=1 

where  the  identity  (a)  from  the  CL  has  been  used.  By  (i)  of  the  Lemma  (section  B.2),  the 
solution  of  the  recursion  is  just  (10). 

Knowing  x k  permits  the  derivation  of  an  important  formula  for  dq-.-q.  which  is  used 
repeatedly  below.  For  unknown  coefficients  cq,  write 

L-k 

dii—ik  =  'y  /  ^  y  ci\---ikJvjt 

t= 0  ji—jtp 

Summing  over  i\  -  ■  -  ik  shows  that  x k  =  Y^t=k  at-k{l)  £q---q  cir  -it-  Comparing  this  with 
Xk  in  (10)  then  gives  at  =  (-1)*,  and  consequently, 

dn-ik  =  E^1)*  E  ^-ikn-jr  (12) 

t- o  ji—jtH 

*i>—  .*jt 
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This  formula,  which  expresses  the  number  of  uniquely  common  names  in  terms  of  the 
number  of  common  names,  is  the  main  tool  in  many  of  the  derivations  which  follow. 

Let  be  the  number  of  exchangeable  names  in  the  dlr..lk  names  uniquely  common 

to  the  state  lists  of  oq , . . .  ,  aik .  Then  the  total  number  of  exchangeable  names  is 

L 

3  =  E  E  jil-ik- 

k — 2  i  i  'i 

Since,  by  (ii)  of  the  PL,  Efo,  .•■**]  =  (1  -  21~k)E[dn...lk\,  then  Eq-i*  E\jh...ik]  =  (1  - 
2 l~k)E[xk\  and  so,  by  (11), 

m  =  m  =  53(1  - 2i~k)  E(- i)t_fe (l)  E  cn-nK ■■■k- 

k= 2  t-k  '  '  *i -•-** 

Reversing  the  order  of  the  t  and  k  summations  finally  gives  the  result  in  (9). 

The  expression  for  /i  can  be  partitioned  to  illustrate  the  marginal  increase  by  adding 
one  extra  agent.  If  /i(L)  temporarily  denotes  the  mean  for  L  agents,  then  (9)  implies  that 
the  marginal  increase  in  the  expected  exchange  rate  of  all  agents  in  going  from  L  —  l  agents 
to  L  agents  through  the  addition  of  agent  is: 

L  _ 

n{L)  —  jj,{L  —  1)  =  A l  ^  ^(— 1)*(1  —  2  )  y  ^  ■  ■  ■  Xit_1-  (13) 

t= 2  *i 


The  expressions  (9)  and  (13)  are  studied  further  in  section  4.5  below. 

It  is  also  possible  to  derive  /i  via  a  similar  argument  but  based  upon  counting  the 
number  of  exchangeable  names  among  those  uniquely  common  to  resource  pools  rather 
than  state  lists.  This  alternative  derivation  is  given  in  section  B.3. 


4.2  A  single  agent 

This  section  solves  the  problem  of  evaluating  m  =  E[Je]  for  the  specific  agent  af.  We 
show  that  the  expected  number  of  exchangeable  names  fie  in  the  state  list  of  agent  at  is 

£/— 1 

He  —  E[Jj]  —  Ai  y  ^(— 1)*+12  t  y  ^  ■■  ■  Xit.  (14) 

t= 1  ii—itjfi 

The  idea  of  the  derivation  of  m  is  the  same  as  for  h  except  that  the  single  index  l  is  frozen 
during  the  process.  If  xtJk  :=  Eij.-ij.,  dth-ik-i  denotes  the  number  of  names  uniquely 
common  to  k  state  lists  which  include  the  list  of  ae ,  we  first  show  that 


(15) 
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and  so,  by  (2), 


EM  =  >~1  E  (-D‘-*+1(tl,)  Ec«.  («) 

t=k—l  '  '  ii  ■■•it 

By  fixing  one  index  as  f  in  (12), 

L-k 

d£ii—ik_1  =  1)  c£ii—ik-ijr-jt 

t=o  jvjty 


so  that  summing  over  ■  ■  -  ik_1  with  the  use  of  (a)  of  the  CL  produces  (15).  Note  that 
since  £fe=.2  xi,k  is  the  number  of  names  in  the  state  list  of  agent  a t  in  common  with  some 
other  state  list  and  therefore  cannot  exceed  the  state  list’s  size,  then 

L  L- 1 

^  n£ 

k—2  t—  1 


always  holds. 

Let  ja1-ik_1  be  the  number  of  exchangeable  names  in  the  names  uniquely 

common  to  the  state  lists  of  a £,  aq,. . .  , Then  the  total  number  of  exchangeable 
names  in  the  state  list  of  a t  is 


L 

~  'y 'y  y  jiii—iic-i- 

k= 2  ii~ik_ijf£ 

Again  by  (ii)  of  the  PL,  E[jhl...ik_ J  =  (1  -  21~k)E[deil...lk_1\  so  that  summing  over  the 
combinations  gives  £ n...4_i  E{jtiv..ik_x}  =  (1  -  2 l~k)E[xt^.  Hence,  by  (16), 

w  =  E[jt\  =  5^(1  -  2i~k)  (l~\ 

k=2  t=k  1 

Reversing  the  order  of  the  t  and  k  summations  gives  the  result  in  (14). 

As  with  all  agents,  the  marginal  increase  in  the  mean  can  be  similarly  deduced  from 
(14).  If  fii(L)  denotes  the  mean  (14)  for  L  agents  present,  then  the  marginal  increase  in 
the  expected  exchange  rate  for  agent  a i  in  going  from  L  —  1  agents  to  L  agents  through 
the  addition  of  agent  ai  is: 


l- 1 

ML) -ML-  1)  =  A,AL^(-l)t+12-f  £  CtLll..,t_, A21  ■•■A1(_1.  (17) 

Again,  both  expressions  (14)  and  (17)  are  studied  in  further  detail  in  section  4.5. 
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4.3  A  syndicate  of  agents 

Consider  now  the  (interpolating)  problem  of  finding  the  expected  number  of  exchange¬ 
able  names  within  the  state  lists  of  a  specified  syndicate  au  =  alx , . .  -  ,  al% ,  2  <  £ >  <  L  -  1, 
of  the  L  agents.  By  re-numbering  the  indices,  it  can  be  assumed  that  ...  ,4  are  the 
first  «  agents.  Let  Is  ,4}  be  these  indices  and  I\IS  be  the  remaining  indices. 

As  in  section  1.1,  let  Ju  be  the  random  variable  of  the  number  of  exchangeable  names 
within  the  state  lists  of  the  agents  aIs  and  let  =  E[JIs]  be  its  mean. 

The  number  ji3  of  exchangeable  names  can  be  written  in  two  parts  as: 
s  s  L—s 

3IS  Jn-rp+EE  S  E  jri—rpii—iiei 


2  r\-rp 

eh 

=: A  +  B 


p= 1  k=  1  ri"’rp  h-"ik 

eh  ei\h 


which  reflects  a  decomposition  into  the  exchangeable  names  inp-cohorts  taken  only  from 
I  (i  e  A)  and  those  in  (k  +  p)- cohorts  containing  at  least  one  agent  from  Is  and  one  agent 
from  I\Is  (i  e  B ).  To  find  p/s,  both  E[A\  and  E[B\  must  be  computed.  Consider  A  first. 
Since  by  (ii)  of  the  PL,  E[A]  =  £^2(1  "  Eri-rP  E[dri...Tp),  the  rewriting  of  the  sum 

L—p 

dry-Tp  =  S  cn-rph-jt, 

t'= 0 

l3T  ii'"  trp 


s—p  L—s 


=  ^^(-i)9+<  E]  £ 

q—0  t=0  h  'lq  h-jt 

els\{ri,--- ,rp}  el\h 

splits  the  summation  over  the  jt>  into  two  disjoint  pieces,  one  within  Is  and  one  within 
I\IS.  Summing  over  the  rp’ s  and  taking  the  expectation  now  gives 

p_2  F0  (=0  ri-rpli-lqji-jt 

s  s  L-s  /\  _  _ 


= i  -  2i_p)  £  E(-1),+<'“p  (!)  E  E 

p=2  q=p  t=0  W  n-rqji-jt 

where  the  CL  and  a  shift  of  p  in  the  ^-summation  have  been  used  Next,  switching  the 
order  of  the  p  and  q  summations  and  using  £^(1  -  ^(-l )p©  =  1  ~  ^  gives 

E[A]  =  EZ('_1)t+<i,(1  ~~  21_g)  £  £  Cri-rqh-itKi  ■  ■  ■  Kq\i  ■  ■  ■  ^ if 


t-0  p= 2 


ri—rq  ii—it 

eh  ei\h 
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Now  consider  B.  As  for  A,  again  write 

L—p—k 

dri—rpii—ik  —  (“l)1  El  cr\---rpii—ikji—jt' 

t'= 0  ji—jt' 

,ik 

—  E  E  (_l)9+t  E  E  cri-rph-ikh-lqji-jt 

9=0  t=0  <i—Z9  ji—jt' 

,rp  ?iiv  >U 

to  split  the  j-summation  over  7S  and  /\/s.  Since  by  (a)  of  the  CL, 


£  £  £  £ 

rv-Tp  h  ■■■lq  h  •••»*  R-dt 
^ri,...,rp  5^1 


cri—rpii-ikh-lqjl-it 


q  +  p\  ft  +  k 
p  )\  k 


£  £ 

n  ii— ijb+t 


-n— r,+p»i— **+< 


substitution  of  the  expectation  of  dn...rpir..ik  into  £?[£?]  gives 


t£(i-21-»-‘)E£Vi)‘-i+»-'’0(‘)  £  £*[*,•■• 

P=1A=1  q=P  t=k  VF/  V  7 


tq*i— id* 


where  the  g  and  i  summations  have  been  shifted  by  p  and  k  respectively.  Interchanging 
the  p  and  q  summations  and  k  and  t  summations  then  produces 


££(-D1+4£(-i)^)£(-D‘(‘)(i 

9=1  t= 1  Lp=  1  k- 1  V  7 


-  2 


l-p-fc'l 


£  £*i- 

n— ii— it 


Ti— r9zi  •••zd 


Evaluating  the  sum  in  the  inner  bracket  as  1  —  2(1  —  2  J)(l  —  2  q)  then  gives 


E[B)  =  E  -  2(*  -  2"t)(1  -  2_<?)]  E  E  C'rx-r^r-i.Ari  ‘  '  ‘  '  '  '  K- 

9=1  t=l  r i  •••r,  ii  ■■■it 


Finally,  combining  the  results  for  E[A\  and  E[B]  produces  the  result 
s  L—s 

n.  =££(-i)t+’(i-21^)  £  £  Cr\--  rqi\  -  it^Ti  '  '  '  ^rq\  1 

9=2  t=0  ri-r, 

eL  6/\/* 

s  L—s 

+  EE(-1)t+9[1“2(1-2"i)(1-2'9)]  E  E  Cri...r,i1...i(Ari---Ar,Ail---Ait 

9—1  Z=1  ri  '"r9  ii  •••if 

6 A  e/\/4 
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or,  equivalently, 


Vis 


— 1)9(1  —  21  q)  CVi  — r,  Ari  ■  ■  •  Ar 


Q= 2 


r  i-r, 

eh 


L-s 


+  EE(-1>'+,2l'‘(1-2^)  E  E  c' 

9=1  <=l 


ri"-rflii—tt^ri  '  '  '  A^Aq  •  ■  ■  Aj(. 


(18) 


ri  ii 

els  ei\h 


Written  in  this  way,  pis  is  the  sum  of  the  exchange  rate  for  the  agents  aIa  considered  in 
isolation  (cf.  9)  plus  a  correction  term  accounting  for  the  interactions  with  the  remaining 
agents  not  indexed  by  I$-  If  the  agents  o jg  form  a  cooperative  syndicate  in  the  sense  that 
they  share  minimal  interest  in  common  resources,  i.e.  CTl...Tp  =  0  for  all  combinations 
ri...r  £  lSj  then  pjs  in  (18)  consists  only  of  the  second  interaction  term. 

As  a  final  check  on  the  above  results,  we  verify  the  means  computed  in  the  previous 
two  sections.  First,  for  the  single  agent  ap  (Is  {^})>  ®  I  ■  ^  this  case  -Ef-A]  must 
be  dropped  from  pt  and  by  setting  q  =  1,  its  only  possible  value  in  E[B],  gives 

L- 1 

W  =  E(-l)t+12_i  E  Ait, 

t—  1  ii—it 

It 

which  is  exactly  the  previous  result  for  a  single  agent  in  (14).  Next,  for  all  L  agents,  set 
s  =  L  (i.e.  Is  —  {1, . . .  ,£}).  In  this  case  E[B]  must  be  dropped  from  p  (no  interaction) 
and  setting  t  =  0  in  E[A]  produces 

L 

p  =  y^(—l)(} (1  —  21  q)  Cri . ..Tq An  •  •  •  Xrq , 

9=2  r  i— r, 

which  is  the  same  as  derived  previously  for  the  system  mean  (9).  Hence,  the  general  result 
for  pis  encapsulates  the  two  previous  cases  (s  =  1  and  s  —  L). 


4.4  The  number  of  common  names 

A  similar  analysis  to  the  previous  section  applies  for  the  number  of  common  names 
to  the  state  lists  of  the  agents  indexed  by  Is.  Let  xIs  be  the  random  variable  of  the 
number  of  names  in  the  state  lists  of  agents  als  common  to  two  or  more  state  lists  and  let 
Xj  :=  E[xj3 ]  be  its  mean,  then 

xi,  =  E  E  *i~*  + EE  E  E  drvrPh  -ik 

p=2  T\  -rp  p=l  k=lri  'rP  ii-  ik 

eh  6L  g i\is 

=:  a  +  b. 

Since  xIs  is  a  count  of  the  number  of  names  appearing  in  the  state  lists  of  agents  aIs  which 
are  common  to  two  or  more  state  lists,  it  tallies  the  number  of  names  in  which  exchange 
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is  possible.  To  find  the  expected  value  xisi  computation  of  E[a]  and  E[b\  is  required. 
A  very  similar  derivation  to  that  of  fijs  above,  but  one  that  does  not  use  the  PL,  shows 
that 

E[a]  =  ^(-l)t+g  (q  -  1)  y!  Cri-rqi1-it^rl  ■■■  ■■■ 

9=2  t— 0  Tl  -Tq 

eis  a\is 

E[b]  =  J2  Y,(-l)t+q  J2  J2  Cr.-rgh-uK  K- 

0=1  t= 1  ri'"r«  *i ---it 

eh  ei\i3 

Consequently,  the  expected  number  of  names  in  the  state  lists  of  agents  ajs  which  are 
common  to  two  or  more  state  lists  is 

Xls  —  ~  1)  r.  Cri-rqXri  •  •  •  Xrq 

9=2  n-rq 

eh 

+  ^2  ^2  ^(_l)f+9  ^2  Cry-Tgii  -itKi  •  •  •  KgXh  '  "  XU  ‘  ^  ^  ) 

9=1  t—  1  r\-  -rq 

€ls  £l\Is 

Again,  this  expression  is  split  as  the  sum  of  the  expected  number  of  common  names  among 
the  agents  ajs  as  if  in  isolation  plus  a  correction  term  to  account  for  the  interactions  with 
the  agents  not  indexed  by  Is. 

One  use  of  the  expression  for  Xh  is  obtain  upper  bounds  on  the  possible  realisations 
of  the  random  variables  J,  Jt  and  Jjs.  This  is  done  by  taking  all  of  the  state  lists  of 
maximum  size  in  which  case  each  is  equal  to  the  corresponding  resource  pool.  For  s  —  L, 
the  term  E[b]  must  be  dropped  to  obtain  an  expression  for  the  expected  number  of  common 
names  among  all  state  lists  (x)-  Setting  t  =  0  in  E[a]  and  every  A ik  =  1  to  give  maximal 
state  list  size  produces 

= B-1)^  -  >>  £ 

9=2  ri--rq 

Xmax  is  the  number  of  names  common  to  two  or  more  resource  pools  and  as  such  is  the 
maximum  possible  number  of  exchangeable  names  among  all  L  state  lists.  In  other  words, 
any  realisation  j  of  J  must  satisfy  j  <  Xmax- 

For  s  —  1,  the  term  E[a]  must  be  dropped  from  the  expression  for  Xi-  Setting  q  =  1 
in  E[b]  and  every  A ik  =  1  to  give  maximum  state  list  size  produces 

=  E(-l)‘+1  E  c*.  -ir 

t=l  ii—itjfl 

Xe, max  is  the  number  of  names  in  the  pool  Pg  which  are  common  to  two  or  more  resource 
pools  and  is  therefore  the  maximum  possible  number  of  exchangeable  names  that  can 
ever  appear  in  the  state  list  of  agent  (if  .  That  is,  any  realisation  jf  of  J(  must  satisfy 

jl  ^  X/,max- 
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For  the  intermediate  values,  2  <  s  <  L  -  1,  and  with  Is  =  {tx,  ■  ,4},  putting  every 
\ik  =  1  in  (19)  produces 

0=2  n-r, 

+  X  zC  X  X]  Cri...rqii-it 

q=  1  t=l  rl'"r9  ii  -it 

eh  <EI\IS 

as  the  number  of  names  in  the  pools  Pjs  which  are  common  to  two  or  more  pools.  Xh,ma.x 
is  therefore  an  upper  bound  for  any  realisation  jj3  of  Jjs  in  that  jis  <  Xh, max  must  always 
hold. 


4.5  Reductions  under  EDCI 

To  easily  visualise  the  behaviour  of  the  mean  with  the  number  of  agents  L,  the  relative 
state  list  size  and  the  number  of  the  selected  agents  aIs ,  it  is  necessary  to  sufficiently 
collapse  the  high  dimensional  dependence  in  the  parameters  and  Xlt.  This  is  done 

using  the  EDCI  assumption  (section  1.1)  which,  by  giving  a  diminishing  common  resource 
interest  with  increasing  numbers  of  agents,  also  allows  closed  form  evaluations  of  the 
multi-index  sums  involved  and  a  subsequent  reduction  to  one-dimensional  polynomials. 

The  EDCI  substitution  Civ..it  =  p*C,  0<p<l,  reduces  the  expressions  for  pIa  in  (18) 
and  xi s  in  (19)  to  simple  (symmetric)  polynomials  in  the  variables  p\t.  For  a  given  index 
set  Is ,  define  the  quantities 


X  =  1-  n  (i-A), 

xs=i-n  (l-pxo, 

iei\h 

leh 

Q. 

1 

1—1 

G 

1 

1- H 

ys  =  i-n(1-A/2), 

l€l\h 

leh 

Zs  =  J2pXi  (1  —  pXm ) 

l€ls  mels\{l} 


First  consider  the  mean  pis.  Substitution  for  Ci r..jt  into  (18)  and  evaluating  the 
resulting  sums  shows  that 


Ms  = 


-  X)(2 Ys  -  Xs)  +  XXS  -  2 (Y  -  X)(YS  -  Xs)  , 


1  <  s  <  L, 


where,  although  derived  only  for  2  <  s  <  L- 1,  the  expression  remains  valid  for  1  <  s  <  L. 
Now  consider  the  expected  number  of  common  names  xi3 ■  Substitution  for  into 

(19)  produces 


Xh 


=  C 


XS  +  (X-  1  )ZS 


1  <s<L. 
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This  expression  is  also  valid  for  the  full  range  of  values  of  s.  Observe  that  the  sizes  of  the 
corresponding  pool  overlaps,  as  computed  in  section  4.4,  are  obtained  by  letting  A/  — >■  1, 
for  all  l,  and  produces 


Xmax  =  [1  -  (1  -  P)L~X  ~  PH  1  -  P)L~X}C , 

Xt,max  =  P[  1  -  (1  ~  P)L~1]C, 

Xla  ,max  =  [1  -  (1  -  PY  -  sp(  1  -  p)L-l)C. 

As  expected,  each  of  these  satisfies  limp_>i  X-,max  =  C  for  the  case  of  a  ‘core’  of  C  common 
resources  and  limp_»o  X-,max  =  0  for  the  case  of  no  common  resources. 

The  efficiency  ejs  Pls/xis  defined  earlier  in  section  1.1  can  now  be  written  as 

(1  -  X){2 Ys  -  Xs)  +  XXS  -  2 {Y  -  X){YS  -  Xs) 
eis  ~  XS  +  (X-  1  )ZS 

To  produce  one-dimensional  plots,  consider  like  agents  for  whom  A i  —  A,  for  all  l,  so  that 

AT  =  1  —  (1  —  pX)L~s,  Xs  =  l-(l-p\y, 

Y  =  1  —  (1  —  p\/2)L-\  Ys  =  1  -  (1  -  pX/2)s, 

zs  =  Sp\{\  -  P\y-\ 

are  all  polynomials  in  pX. 

For  s  =  L,  II  =  {1, . . . ,  L}  designates  all  agents  and  the  values  X  —  Y  =  0  should  be 
used  in  the  above  expressions.  Hence,  p  =  (2 Yl  —  Xl)C  and  x  =  {Xl  ~  Zi)C.  If,  as  in 
section  4.1,  p(L  —  1)  denotes  the  mean  number  of  exchangeable  names  for  L  —  1  agents  in 
isolation,  then  the  relative  marginal  increase  by  adding  agent  ai  is 

A ip(L  -  1)  _  x  Xl- i  -  Yl- \ 
p(L-l)  ~  P  L  2YL-\  -  Xl-i  ' 

For  the  further  reduction  A i  —  A,  for  all  /,  Figure  4(a)  plots  the  mean  p/C  while  Figure 
4(b)  plots  the  relative  margin  A lp{L  —  l)/p(L  —  1),  both  as  functions  of  pX  for  various 
values  of  L.  Whereas  the  mean  increases,  the  relative  margin  decreases  as  a  function  of 
both  pX  and  L.  Figure  5(a)  plots  the  efficiency  e  =  (2Yz,  —  Xl)/{Xl  —  Zi)  for  various 
values  of  L  as  a  function  of  pX.  Figure  5(b)  plots  efficiency  contours  in  pX-L  space.  These 
are  generated  by,  for  each  efficiency  level  e,  solving  e(pA)  =  e  for  pX  for  each  fixed  L.  Note 
that  certain  high  levels  of  efficiency  are  not  attained  without  sufficiently  many  agents  as 
is  indicated  by  the  failure  of  the  corresponding  curves  in  Figure  5(a)  to  ever  reach  those 
efficiency  values. 

For  s  —  1,  lx  =  {£}  designates  the  single  agent  ae  and  Z\  —  X\  =  2Yi  =  pX^.  In  this 
case,  the  mean  is  pe  =  pX(YC  and  the  expected  number  of  common  names  is  xe  —  pX^XC. 
If,  as  in  section  4.2,  pe(L  —  1)  denotes  the  mean  for  agent  a £  among  L  —  1  agents  in 
isolation,  then  the  relative  marginal  increase  through  adding  one  more  agent  a,L  is 

Alpi(L-I)  ,  x  l-Y(L-l) 

ML- 1)  . *P  L  Y(L-  1)  ’ 
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Figure  f:  (a)  A  plot  of  p/C  for  L  =  2, ...  ,10  (bottom  to  top);  (b)  plot  of  the  relative 
margin  A Lp(L  -  1  )/p(L  -  1)  versus  p\  for  L  =  3, . . . ,  11  (top  to  bottom). 


Figure  5:  (a)  A  plot  of  the  efficiency  e(p\)  for  L  =  3, . . .  ,10  (bottom  to  top);  (b)  efficiency 
contours  in  L-pX  space  for  e  —  0.95,0.90,...  ,0.55  (top  to  bottom). 


where  Y{L  -  1)  denotes  Y  defined  only  over  the  agents  a i,...  ,aL- 1-  Figure  6(a)  plots 
the  ratio  pe/(pX(C)  and  Figure  6(b)  plots  the  relative  margin  ALpe{L  -  l)/pe(L  -  1),  for 
various  values  of  L,  as  functions  of  pX  under  the  further  reduction  Xm  =  A,  for  all  m  ^  l. 
Again  the  mean  increases  and  the  relative  margin  decreases  as  functions  of  both  pX  and 
L.  The  efficiency  e£  =  Yj X  is  very  similar  in  shape  to  e  as  a  function  of  both  pX  and  L. 
However,  ee  >  e,  and  if  Am  =  X(  =  X  for  all  m,  then  the  peak  difference  remains  bounded 
according  to  0.02  <  maxL  maxpA(e£  -  e)  <  0.06  as  illustrated  in  Figure  7. 

Next  consider  the  convergence  of  pjs  to  p.  and  e/s  to  e  as  s  increases  from  1  to  L. 
Again  by  putting  A;  =  A,  for  all  /,  Figure  8  shows  plots  of  the  mean  pIs  in  (a)  and  the 
efficiency  e/s  in  (b)  for  L  =  10  and  s  =  1, . . .  ,  10.  Observe  that  while  pIs  increases  from 
the  individual  mean  p£  to  the  system  mean  p  with  increasing  s,  er,  decreases  with  s  from 
the  individual  efficiency  et  to  its  system  value  e  for  all  L  agents.  Hence,  when  considered 
as  a  group,  any  syndicate  of  s  agents  enjoys  a  higher  efficiency  but  a  lower  mean  than  the 
whole  collection  with  monotone  convergence  to  the  A-agent  system  values  as  s  increases. 
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Figure  6:  (a)  A  plot  of  [j.e/(p\eC)  (bottom  to  top);  ( b )  the  single  agent  relative  margin 
A lPi{L  —  1  )/pe(L  —  1)  as  functions  of  pX  for  L  —  3, . . .  ,10  (top  to  bottom). 


Figure  7:  A  plot  of  the  efficiency  difference  ee  —  e  as  a  function  of  pX  for  the  values 
L  =  3,...  ,10. 

Figures  4,  5,  7  and  8  are  indicative  performance  curves  for  networks  of  like  agents.  Of 
course,  if  the  many  parameters  are  known,  the  exact  formulae  from  section  4  can  be  used, 
however,  the  simple  one-dimensional  plots  here  provide  intuition  as  to  how  performance 
should  change  with  both  the  number  of  agents  present  and  the  fractions  of  their  pools 
represented  within  state  lists. 
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Figure  8:  Plots  of  the  means  Hh  in  (a)  and  the  efficiencies  eIs  in  (b)  as  functions  of  pX 
for  the  values  s  —  1, . . .  ,10  with  L  =  10. 


5  Conclusion 

This  report  has  developed  a  combinatorial  probability  model  for  the  expected  exchange 
rate  of  resources  when  L  >  2  unbiased  agents  meet  in  a  crossing  network.  The  unbiased 
assumption  has  allowed  the  derivation  of  explicit  formulae  for,  in  simple  cases  (L  =  2,  3), 
the  exchange  distribution  and  certain  higher  moments  of  the  distribution.  In  particular, 
for  T  —  2  this  distribution  is  favourably  (positively)  skewed,  but  only  significantly  so  for 
small  values  of  the  mean  and  variance,  a  property  which  rapidly  diminishes  as  state  list 
size  increases. 

For  general  L,  the  mean  of  this  distribution  (i.e.  the  exchange  rate),  has  been  found 
for  all  of  the  agents,  a  single  agent  and  an  arbitrary  syndicate  of  agents.  Through  a 
suitable  reduction  of  the  number  of  defining  parameters  (edci),  reasonable  intuition  can  be 
gained  for  the  behaviour  of  crossing  networks  through  the  development  of  one-dimensional 
expected  outcome  plots  as  were  detailed  in  section  4.5.  For  example,  while  the  means  p 
and  pi  are  increasing  functions  of  L  and  state  list  size,  the  relative  margins  A Lp/p  and 
A LPe/Pi  are  both  decreasing  functions.  The  individual  efficiency  et  exceeds  the  system 
efficiency  e,  however,  both  are  increasing  functions  of  relative  state  list  size.  Such  intuition 
(see  Figures  4  -  8)  is  useful  for  administrators,  regulators  and  users  of  crossing  networks 
for  both  the  prior  prediction  (the  participation  decision)  and  post  evaluation  (the  re-use 
decision)  of  exchange  outcomes. 

For  the  general  case,  much  remains  to  be  done.  Having  here  only  derived  the  mean 
exchange  rate,  the  question  of  higher  moments  remains  open.  For  example,  computation 
of  the  variance,  useful  in  identifying  outlying  events,  requires  an  appropriate  method  to 
compute  the  covariances  of  the  common  name  variables  Cjr..Zf  defined  in  section  1.1.  As 
was  seen  with  the  case  L  =  3,  this  seems  not  to  be  a  simple  task  if  first  finding  the  joint 
distributions  is  required.  This  hints  that  perhaps  an  alternative  approach  to  the  derivation 
of  the  second  and  higher  moments  is  required. 

A  second  major  open  problem  is  the  dropping  of  the  unbiased  assumption  itself.  The 
uniform  prior  distribution  on  both  appearance  of  resource  names  in  state  lists  and  the 
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require/release  intention  for  each  resource  is  the  ‘easy’  case  to  consider  and  may  rarely 
hold  in  reality.  A  model  therefore  needs  to  be  developed  for  given  (non-uniform)  discrete 
distributions  for  these  quantities.  While  in  this  case,  simulation  provides  an  alternative, 
the  dimensionality  of  the  parameter  space  describing  L  agents  is  likely  to  be  too  large 
to  gain  a  simple  picture  of  the  relevant  moments  of  the  exchange  variables  as  has  been 
gained  here  in  the  unbiased  case.  An  analytical  result  with  simple  reductions,  as  was  done 
in  section  4.5,  would  thus  be  more  useful. 

A  number  of  combinatorial  identities  (cf.  sections  3.2.2  and  A.l)  have  resulted  from 
the  analysis  of  the  exchange  model  -  particularly  from  the  case  L  =  3.  It  remains  to  be 
seen  if  these  are  in  fact  new  or  can  be  recast  as  known  identities,  perhaps  as  generalised 
convolutions,  and  therefore  fit  into  an  existing  framework  [6].  This  remains  a  topic  for 
further  investigation. 
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Appendix  A:  Some  Computational  Details 


This  appendix  contains  quite  a  bit  of  gory  detail.  It  provides  the  more  computationally 
intensive  part  of  section  3.2.1  (L  =  3)  by  deriving  an  expression  for  T;m(r )  -  this  result 
appears  in  section  A. 2.  The  results  and  formulae  derived  here  are  more  general  than  are 
required  for  this  task  but  are  of  interest  in  their  own  right  as  giving  methods  of  computation 
and  identities  for  certain  sums  in  terms  of  the  DG-distribution  P2(C;  r)  discussed  in  section 
2.2.  Once  these  general  results  are  established,  it  is  an  easy  matter  to  specialise  them  to 
recover  the  value  of  T/m(r). 


A.l  General  identities 


We  begin  by  considering  sums  which  involve  the  same  terms  as  the  DG-distribution 
p2  (C;  r)  but  summed  against  powers  of  the  index.  Simple  versions  of  these  sums  arise  in 
taking  the  expectation  E\cimciq }  of  common  name  variables  in  section  3.2.1. 


Proposition  1.  For  ni  <  Ni,  C  <•  B ,  A  <  B  and,  s  —  0,1,2,...,  the  identities 

Kit?) ©(?--;) 

Ss{r):=2^k  ~ 


ffi  (?) 


=  £ 

71  =  0 

5 


(-1)"  Mr 


cn,p(A)  (~Ar-  •  r)PP2(C;  71/  -  n,  A;  Nt  -  n,  B;  r) 


£  !§ik  £  C  W-M  V2 (C;  n,  -  n,  A  -  n;  N,  -  n,  B; r) 


,=0  nl^n  ^ 


(Al) 

(A2) 


hold,  where 


P2(C;nhA-,NhB-,r)  :=  ^ 

k 


©  (?) 


is  the  DG- distribution,  Ar  is  the  forward  difference  operator  in  r  and  cn>p(A)  and  aq  are 
expansion  coefficients. 


Before  proceeding  with  the  proof,  two  points  should  be  noted. 

Remark.  The  degenerate  case  B  =  A  =  Cim  (and  s  —  1)  is  the  one  ot  interest  for  the 
calculation  of  T,m(r)  in  section  3.2.1.  When  B  =  A,  the  DG-distribution  reduces  to  the 
hypergeometric  distribution: 


(C)  (N'-Cj 

P2(C;nl,A-,Nl,A;r)  =  HriC-^Nt)  :=  W  M  (A3) 


for  which  the  results  (Al)  and  (A2)  remain  valid. 
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Remark.  By  defining  the  nonnegative  quantities  pk  :=  (^)  (nj)  x  (r)  (c-r)/(c)>  so 

that  P2(C;r)  =  the  defining  sum  for  Ss{r)  in  the  proposition  can  be  regarded  as  a 

multiple  of  the  s-th  ordinary  moment  of  the  discrete  (in  k )  distribution  pij P2{C',r). 


Proof  of  Proposition  1.  The  required  results  are  derived  by  a  straight  forward  application 
of  the  “snake  oil”  method  [7].  Defining  the  generating  function  Ss(z)  :=  Ss(r)zr ,  then 
since  the  hypergeometric  terms  (^)  (clr)/(c)  are  Senerated  by  2F\{—C,  — k ;  —  B\  1  —  z), 


n  /  \  Sf'  i  s  (fc)  (nj-A:)  ( ~C)t(-k)t  (1  ~  zf 

,w  "V  o  v  <-*>«  <! 


V'  ;  J[/ .]  (*:)  (nj-fc)  (“CM-1)*  (1  ~  ZY 

=  w[  ]t  O  <-*><  « 

V-r, ,  (i)  (n!-k)  (-C)t(-iy  (1  -  *)* 

=  2^2^anMt)2^W*+”--(Nfj - (Z ~B)t - 


t  n— 0 
s 


t\ 


CfR-n  [nd*+n  (-C)t(-iy  (1  -  ZY 

:)[1,+"w)<«  (~b)‘  — 

_  [ nl]n  Y^  /.\r  a  .1  -n]t(-A)t  (1  ~  z) 

~  2^  fA/,1  zA°n>s^^  t\ 


n=Q  t 


(-l)nMr 


£  n,M)» 


EE  Cn,p(Z)tP 


[AT/-n]t(-B)<  t\ 
(~C)t[ni  -  n]t{-A)t  (1  -  z)1 


t  p= o 


[Ni  -  n\t(—B)t 


t\ 


(— l)n[n(]n  t  w(  <  d  \p  (~C)tin  ni)t{  A)t{l  z)1 

=  E^pjrE^)((--1)*J  E-  {n-mlzsu — “ 


-  E  ^T1  E  ((.  -  !>£)'.*(-* »  -  -*»  -  N„  -B;  1  -  ,). 

n=0  L  l>n  p= 0 


d  \P 


Since  the  operator  (z  —  l)d/dz  on  generating  functions  corresponds  to  the  operator  —A r-r 
on  coefficients,  and  ^F2(—C,—m,—n;-M,--N-,l  —  z)  —  Yhr  P2{C;rn,n;  M,  N\r)zr ,  the 
last  identity  produces  the  required  result.  A  similar  chain  of  identities  but  with  the  split 
[A]f+n  —  [A  —  n]t[A]„  produces  the  equality 

=  E  [j$t  £  O' "•  (<2 "  ^s) ' : A(~c' ~n‘ +  -A  + w' ~B' ~N‘ + m 1  _  ■ 2) 


which  is  the  second  form  of  the  sum  in  (A2). 

To  make  use  of  either  of  the  forms  (Al)  or  (A2)  requires  expressions  for  the  coefficients 
Cn,p{A)  and  aq.  The  terms  an,s{t)  are  defined  through  ks[k]t  —  YHn=oan,s{t)[k]t+m  while 
the  terms  (— l)nc„)P(A)/n!  are  the  coefficients  of  tp  in  a7ltS(t)[A  —  t]n.  The  degree  s  —  n 
polynomials  antS{t)  satisfying  the  recursion 

ao,s(t)  =  t  a0,s-i(t),  n  =  0, 

®7i,s(f)  =  ®n—  l,s  —  1  "b  1)  T  t  0,n,s  —  l(t),  7T.  >  0. 
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Consequently,  aotS(t)  =  ts,  ai^s(t)  —  (l+t)s  ts,  a2,sW  —  £g=o  tq[{2+t.)  q  (1+^)  9  ] 

and  generally, 

s—n 

an,s(t )  =  i9an-l,s-l-g(l  +  <)> 

9=0 

from  which  there  follows  the  closed  form  expression 


“»..(<) = ^  d-^C)  (t + «>* = s  £  “ « where 

9=0  ?=n 

Here,  the  introduction  of  the  coefficients  aq  is  achieved  by  expanding  the  (t  +  q)s  term. 
To  finally  obtain  c„;P,  use  the  explicit  form  of  anyS(t)  to  expand  a„)S(£)[A  -  t]n  to  get 

v 

Cn,p  =  Pm'ip-mi  p  =  0,  . . .  ,S, 
m— 0 

A»  =  E  [i  m  = 

l—m  L  J 

7 9  =  as-9  f  ^  >  g  =  0, . . .  ,  s  —  n, 


with  aq  as  above  and  [”]  the  first  kind  Stirling  number  [7].  Observe  that  the  two  forms 
(Al)  and  (A2)  originate  from  the  alternative  splits  [. A]l+n  =  [A\n[A-n]t  =  [4]t[A-i]n.  □ 

For  reference,  we  compute  and  list  a  few  of  the  coefficients  cntP(A)  as  (s  +  1)  x  (s  +  1) 
matrices  with  row  index  n  and  column  index  p,  both  beginning  at  zero  and  ending  at  s. 
Of  course,  for  s  —  0,  the  only  nonzero  coefficient  is  eo,o(4)  =  1-  Below  are  the  matrices  of 
coefficients  for  s  =  1, . . .  ,  4: 


s  =  1  : 


s  =  2  :  — 


s  =  3  : 


-A 

2A(A-1)  -2 

0 

- A 

6 A  (A-  1) 


-(2A-1)  2 

2(2A-1)  2 


-(34-1) 

6  (A2  -3A  +  1) 


-3  (4-1) 
—12  A 


-6  4  (4-1)  (4-2)  6  (3 A2  -6A  +  2)  -18(4-1)  6 


s  =  4  :  — 


14  4  (4- 1) 


-(4  4-1)  -2(3  4-2)  -2(2  4-3)  22 

2  (l2  42  — 26  4+7)  2  (6  42  — 30  4+19)  -12(2  4  -3)  22  3 


—36  4  (4  —  1)  (4  -  2)  -12(24-3)  (42  -64+2)  12  (6  42  -21  4+13)  -36(24  -  3)  23  3 
24  [4]4  -48(24-3)  (42-34+l)  24  (6  42-18  4+ll)  -48(24  -  3)  23  3 
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Whereas  Ss(r)  gives  the  ordinary  moments  oi pk/ P2{C',r)  (see  the  Remark  above),  the 
factorial  moments  can  also  be  considered  by  defining 

*r  :=?  W- Cl  (IT' 


The  defining  relations 


=  and  *‘  =  E{i}W' 

i=0  L*J  i= 0  k  } 

for  the  first  and  second  kind  Stirling  numbers  [7]  imply  the  corresponding  identities 

cs(r) =  i Si^  and  Ss^  =  S*!  {*}  Ci^- 
2=0  <■  -*  2=0 

These  in  turn  lead  to  a  companion  result  to  (Al);  namely, 

a s(r)  =  ”T  wT  S  (~&r  ■  r)PP2(C;  nt  -  n,  A;  Nt  -  n,  B-,  r), 

where  the  coefficients  (— l)ncnp(A)/n!  are  now  the  coefficients  of  tp  in  aHtS (t) [A  —  t]n  with 


~  /  \  1  v — >  S  ,  v  1  \  ^  ® 

072,4*)  ==  2^  r  a"’r^  =  1  Jt 

‘  T=72  L  J  <7=22  \ 


and  the  same  a„>r(t)  as  above.  The  first  few  sets  of  coefficients  c„jP(A)  are,  for  s  —  1, . . .  ,  4: 


s  —  1  : 


8  =  2  : 


|(  o  ~2A  2) 

\2  (A  -  1)  A  -2  (2  A  -  1)  2/ 

/  0  2  -31 

1  0  3  A  — 3  (A  + 1)  3 

3!  0  6  (A-l)  A  -6  (2 A-l)  2-3 

\-6  (A -2)  (A-l)  A  6  (3  A2  -  6  A  +  2)  -18  (A  -  1)  2-3 

/  0  -6  11  —6  1  ^ 
0  -8  4  4(3  4+2)  -4(4+3)  22 

—  0  -12(4-1)4  12  (42+4  — l)  -244  22  3 

4'  0  -24  (4-2)  (4-1)  4  24  (3  42  — 6 4+2)  -72(4-1)  23  3 

\24[4]4  -48(24-3)  (42-3  4+1)  24  (6  42-18  4+ll)  -48(24-3)  23  3/ 


A. 1.1  Recursive  evaluation  of  (Al) 

One  alternative  for  the  evaluation  of  the  RHS  of  (Al)  is  to  use  the  identity  (for  integer 
n  >  0) 

Hr(A,m  -  n.JV,  -n)  =  ~^-Hr(A,nl,Nl)  (A4) 

\nl)n\^l  A)  n 
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directly  in  the  hypergeometric  term  in  P2{C\ni  -  n,A;N{  -  n,B-,r).  By  then  using  the 
expansion  [k-l]n  =  (-l)n  E"=0(“l)<<Tn-<fc'>  where ai  :=  0+n~ 3.),  0  <  *  <  n, 

are  the  elementary  symmetric  functions  in  l,  l  +  1,  ■■*,/  +  n  —  1,  produces  the  equality 


P2(C;  n/  -  n,  >1;  IVj  -  n,  5;  r)  = 
(iV;  -  A)n 


M)n 

{n0n{Nl  -4)n 


£(-i  )Vi5iW- 

z=l 


(A5) 


When  used  in  (Al),  (A5)  expresses  5s(r)  in  terms  of  S'j(r),  1  <  i  <  s,  and  so  leads  to  a 
five  step  procedure  for  the  recursive  evaluation  of  Ss(r): 


1.  Use  (A5)  to  write  P2(C;rp  -  n,A;Ni  -  n,B;r)  in  terms  of  P2{C-,ni,  A;  Nt,  B;r)  and 
Sz(r),  l<i<s; 

2.  Substitute  for  P2(C;n,  -  n,A;Nt  -  n,B;r)  in  (Al)  of  the  Proposition  to  obtain 
an  order  s  difference  equation  for  Ss{r )  involving  P2(C;ni,  A-,  Nt,  B-,r)  and  Sj(r), 
1  <  i  <  s  —  1; 

3.  Using  the  correspondence  (1  -  z)d/dz  A r  ■  r  between  generating  functions  and 

coefficients,  convert  the  difference  equation  to  an  order  s  ode  for  the  generating 
function  Ss(z)\ 

4.  Solve  this  ode  for  Ss(z)  modulo  C(3F2(1  -  z))  =  0,  where  C  is  the  hypergeometric 
ODO  for  the  generating  function  of  P2(C;  ni,  A\  Nt,  B;  r); 

5.  Pass  from  the  generating  function  >Ss(.z)  to  obtain  the  coefficients  Ss(r )  in  terms  of 
the  Sj(r),  1  <  j  <  a  -  1,  and  P2{C;nh  A;  NhB;r). 

To  illustrate,  these  steps  are  implemented  for  s  =  1.  Performing  steps  1  and  2  gives  the 
first  order  difference  equation 

Ar{rPi(r)}  +  NiS\(r)  =  (A  +  nt  -  Ni)Ar{rP2{r)}  +  ri[AP2{r), 

where  the  abbreviation  P2(r)  =  P2(C;  nh  A;  Nt,  P;  r)  has  been  made.  Step  3  then  gives  the 
first  order  ODE 

(1  -  z)S\{z)  +  NiSi(z)  =  {A  +  nt  -  Ni){l  -  z)P^(z)  +  niAP2{z). 

Applying  an  integrating  factor,  this  is  the  same  as 

{(1  -  z)~Nl Si{z)}'  =  {(1  -  z)-N,[aP2(z)  +  (Po  +  f3iz)P2{z)  +  yz{l  -  z)P%{z)]}' 
-7(1  -z)-n'~1C{P2{z)), 


where 


a  - 


niA 


Ni- 

00  =  ~  1  + 


C’ 

A  +  ni  —  B  —  1 
Ni-C  : 


-1 


Pi  —  1  + 


1  -  A-ni 

Ni-C 
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and  C  is  the  ordinary  differential  operator  given  by 

C{u)  =  z(  1  -  zfu'"  -  (1  -  z)[  1  -Ni-B-{3-C-ni-A)(l-  z)]u" 

+  [N[B  —  (1—  C  —  fii  —  A  -C  niC  +  niA  +  AC)(1  —  z)]u  —  niCAu. 

Since  the  generating  function  P2{z)  =  3F2(—C\  —ri[,  —A;  —Ni,  —B\  1  —  z)  lies  in  the  kernel 
of  C  [5],  then  integrating  gives  the  relation  (step  4) 

Sx{z)  =  aP2(z)  +  (A,  +  Piz)P^(z)  +  yz(l  -  z)P2{z). 

Finally,  returning  to  the  generating  function’s  coefficients  (step  5)  gives  the  relation 

-  [r+(n,A-C~r)]^r)  "  <r  +  ^  +  (A6) 

This  expresses  S\{r)  as  a  simple  quadratic  (in  r)  combination  of  P2{C\  r)  with  index  shifts 
of  0  and  1. 


A. 2  Evaluation  of  7}m(r) 

This  section  gives  two  derivations  of  the  value  of  Tim(r)  corresponding  to  the  two 
techniques  above  for  the  evaluation  of  Si(r).  The  case  of  interest  is  B  —  A  and  hence,  as 
noted  in  (A3),  the  reduction  P2{C\  A,ni;  Ni,  A;r)  =  Hr(C,ni,  Ni)  holds. 


A. 2.1  Method  1 


For  s  =  l,  the  coefficients  co,o  =  0,  co,i  =  1,  c^o  =  —A  and  c\t i  —  1  from  above  can  be 
substituted  into  (Al)  to  give 

5i(r)  -  -Ar  ■rffr(C,n/,JV,)  +  (-l)|{— A+  (-Ar -r)}Hr(C,ni  -  l,Nt  -  1).  (A7) 

Now  the  hypergeometric  probabilities  Hr  satisfy  (put  n  =  1  in  (A4)) 

^Hr(C,nt  -l,Nl-l)  =  ^Fr(C,ni,A) 


and 


-A  rrHr(C,  n,  -  k,  Nt  -  k)  =  ^  +  r  -  IT'  +i~Hr(C>  nl  -k,Nt-k), 


Ni  -  C  +  r  -  nt  + 

for  k  =  0, 1, ... ,  which,  with  k  —  1  simplifies  (A7)  to 


S\(r)  =  - —  + 

Setting  A  =  Cim  finally  gives  the  required  result 


}Hr(CtnhNi). 


Ttm{r )  = 


ni  —  r  rNi  —  niC 
Nt-C  +  Nt-C 

:d  result 

Cirn  {ri[  -r)  +  rN[  -  ntC  £)  (n‘,-r) 


Nt-C 


o 
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A. 2. 2  Method  2 


From  the  expression  (A6)  for  Si(r)  derived  by  the  ODE  method,  we  need  only  set 
B  =  A  =  Cim  and  substitute  Hr  for  P2(r).  Since  Hr  satisfies  the  recursion 

^c,n„Nl) = (r+1)((cw;_r)c(n_;;,lTT)  wm 

then  substitution  into  the  RHS  of  (A6)  produces  the  same  result  for  T/m(r). 
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Appendix  B:  Derivations  and  Identities 


After  deriving  the  three-agent  common  name  distribution  in  section  B.l,  section  B.2  of 
this  appendix  lists  some  elementary  but  useful  combinatorial  identities,  including  what  is 
called  the  Combinatorial  Lemma  (CL).  Finally,  in  section  B.3,  an  alternative  derivation 
of  the  mean  exchange  rate  for  all  agents  is  given. 


B.l  Derivation  of  P(c/m,c/g)  for  L  =  3 

In  this  section  the  LTP  and  Bayes’  rule  are  applied  to  derive  the  joint  distribution 
P(cim,  ciq )  of  common  names  as  stated  in  section  3.2.1.  Define  the  notation  P(k  of  Cim  in  Z) 
to  be  the  probability  that  k  of  the  C{m  names  common  to  pools  Pi  and  Pm  are  in  the  state 
list  of  agent  o/,  etc.  By  the  LTP, 

P(clm\cig)  =  Y^P(cim\k  of  cim  in  l)P(k  of  Cim  in  l\C[q), 

k 

where  the  first  probability  in  this  sum  is  (cfm)  (n^-c^/Cl)'  Considering  the  second 
term  in  the  sum  and  applying  Bayes’  rule  gives 

P(k  of  Cim  in  l\clq  Ciq  in  Z  and  q)  — 

P(k  of  Cim  in  Z) 

P(q,  of  C„  in  (  and  #  of  C,m  in  /)  x  p(ci?  of  Qq  ,  and  (B1) 

where  the  ratio  of  the  two  terms  on  the  RHS  is  just: 

P(k  of  Ctm  in  l )  _  (6fcm )  (7Vn,-A:m) 

P{ciq  of  Ciq  in  l  and  q)  (^)P2(C'(9;  ctq) 

Now  consider  the  first  term  on  the  RHS  of  (Bl)  and  again  apply  the  LTP  (twice)  to  get 
P(ciq  of  Ciq  in  l  and  q\k  of  Cim  in  l )  = 

^  P(ciq  of  Ciq  in  l  and  q\r  of  C  in  l)P(r  of  C  in  l\k  of  C'irn  in  l) 

r 

=  ^  P(ciq  of  Ciq  in  l  and  q\s  of  Ciq  in  l)P{s  of  Ciq  in  l\r  of  C  in  l ) 

r,s 

x  P(r  of  C  in  l\k  of  C[m  in  l). 

The  first  term  on  the  RHS  is  ( s  )  (  Nq~s)/(^q)  while  applying  Bayes’  rule  to  the  second 

' Ciq /  \Tlq  C\q  /  '  \Tlq  / 

gives 


P(s  of  Ciq  in  Z|r  of  C  in  /)  = 


P(s  of  Ciq  in  l) 
P(r  of  C  in  Z) 


P(r  of  C  in  Z|s  of  Ciq  in 


Z). 


(B2) 
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Now  the  numerator  of  the  ratio  in  the  RHS  of  (B2)  is  (C^)  ©-©/©  while  the  denom¬ 
inator  is  (£) inllr) / in1,) ■  The  remaining  terms  in  the  product  are  simply 

fS) 

P(r  of  C  in  l\s  of  Ctq  in  l )  =  r 

(  c) 

and 

(Clm 

P(r  of  C  in  l\k  of  Clm  in  l)  =  ■ 

V  c) 

Finally,  putting  all  of  these  terms  together  and  multiplying  through  by  the  probability 
density  P{clq  of  Ciq  in  l  and  q)  =  P2(Ciq;ciq)  gives  the  desired  expression  for  P{clm,clq) 
in  (5). 


B.2  Elementary  identities 

The  derivation  of  the  means  in  section  4  relies  on  simple  identities  involving  multi- 
indexed  summations.  Although  only  ever  using  part  (a)  in  section  4,  the  results  below 
are  listed  for  reference.  Parts  (b)  and  (c)  have  applications  in  the  following  section  and  in 
reductions  of  the  means  as  in  section  4.5. 

(Joj^iLinatorial  Lemma  (CL).  Fov  any  multi-indexed  quantities  —  t  L ,  the 

following  identities  hold: 

(a)  For  0  <  r  +  k  <  L, 

V.  T.  Qil-ikh-Jr  =  1  X/  Qh-ik+r'l 

il—ik  jl—jryf  *1— »fc+r 

ill---  ;ifc 

(b)  If  s(\[k],r)  is  the  r-th  symmetric  polynomial  in  variables  A^,...  ,  \k ,  then  for  r  < 
k  <  L  and  0  <  t  <  L  —  k, 

(k  J  ~t>  _  ^ 

n—ifc  n—ik+t 

*i>—  ,4 

(c)  If  s(\[k],t)  is  the  t-th  symmetric  polynomial  in  variables  A21,...  ,  A  *fc,  then 

s(A[A:])^)  =  t  =  0,---,k. 

The  proofs  of  these  facts  follow  by  simple  counting  arguments  on  the  possible  number 
of  combinations.  The  following  elementary  result  is  also  noted. 
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Lemma.  For  the  finite  sequence  po, . . .  ,pn,  then  the  identities 

k-r  V  7 


hold.  Additionally,  for  any  positive  integer  s  and  with  b *  defined  as  above,  then 

ksbk  =  ([)  (“Ar  •  r)spr 

r- 0  ^  7 

(-&r-ryPr=£(-irfik)k’bt 

k-r  '  7 

where  Ar  is  the  forward  difference  operator. 

This  Lemma  follows  from  an  application  of  the  ‘snake  oil’  method  [7].  When  the 
sequence  { pT }  is  a  probability  density,  then  k\bk  =  are  the  factorial  moments  and 
hence  pr  is  written  in  terms  of  these  moments.  This  gives  another  easy  way  of  recognising 
the  factorial  moments  of  the  DG-distribution  from  the  expression  for  Pt(C]c)  in  section 
2.2. 


B.3  Alternative  derivation  for  all  agents 

An  alternative  derivation  for  the  expression  for  p  in  section  4.1  which  is  based  upon 
a  consideration  of  the  overlaps  of  resource  pools  rather  than  state  lists  is  detailed  here. 
Although  this  approach  is  a  little  more  involved,  it  is  included  as  a  matter  of  side  interest 
since,  historically,  it  provided  the  first  derivation.  Some  slight  notational  differences  to 
the  main  body  of  the  paper  are  used. 

Let  ji  j  ...jf  denote  the  number  of  exchangeable  names  in  the  state  lists  of  agents  an , . . .  ,  ap 
among  those  names  uniquely  common  to  the  pools  Pn, . . .  ,Pit.  Then  the  total  number  of 
exchangeable  names  is 


i  =  (B3) 

t=2  i\ ■■■it 

The  idea  here  then  is  to  add  up  the  number  of  exchangeable  names  in  the  uniquely  common 
names  to  all  pools  rather  than  state  lists.  By  an  argument  analogous  to  the  derivation  of 
(12),  the  number  Cq...^  of  names  uniquely  common  to  pools  Pq , . . .  ,  Pi,,  is  just 

CLit=E(- 1)‘  E  (B4> 

t= 0  jl— 3 1 p 

ii,...  ,ik 
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which  also  gives  the  upper  bound  on  jix-ik ■  From  this  expression,  for  each  fixed  l  = 
1, . . .  ,  L,  if  0^  denotes  the  number  of  names  in  pool  P(  common  to  all  other  pools,  then 
since  the  total  overlap  of  pool  P(_  with  all  other  pools  cannot  exceed  its  size,  the  inequalities 

c?  :=  E  E  <4U  =  E(-Dl+1  E  £  N‘ 

fc— 1  i\--ik  k—1  ii ••• ik 

must  hold  for  each  £  =  1, . . .  ,  L.  This  provides  a  consistency  check  on  the  specification  of 
the  overlap  parameters  Civ..ik. 

To  derive  p  =  E[j],  first  requires  E[jir..i J.  To  this  end,  we  first  show  that 

EUi, --ij  =  El-1)^1  -  2l"‘MAW,<)  <B5) 

t= 2 

where  s(A[£],  t)  is  the  f-th  symmetric  polynomial  in  A q ,  •  •  •  ,  A?t; .  To  minimise  the  subscript 
notation,  it  suffices  to  prove  this  for  k  -  L  (which  can  be  replaced  by  any  smaller  subset 
of  the  L  agents).  Among  the  C[[..L  =  C  names  uniquely  common  to  pools  Pi,...  ,Pl,  let 
c'  be  the  number  of  these  names  common  to,  and  d!iv..ik  the  number  of  these  names 
uniquely  common  to,  the  state  lists  of  agents  ah,...  ,aik,  and  let  j'iv..ik  be  the  number 
of  exchangeable  names  among  the  d'h...ik  uniquely  common  ones.  Then  by  (ii)  of  the  PL, 
E\j’h..,k\  =  (1-2 l-k)E[d'n...ik}.  Now,  since 

d'  =  c  ■ 

U---U  *1  •••** 

then,  through  an  argument  the  same  as  in  section  4,  and  using  (c)  of  the  CL,  x'k  := 
En  ..4k  EK  -J  is  g^en  by 

xi  =  cE(-1)r'‘(l)*<AW'r)- 

r=k  '  ' 


L-k 

-E  E  ^ h—ikji—jt 

t=l  ji-jt? 

*1  !■■■  >4 


Now,  writing  the  decomposition 

L 

jl—L  —  ^  ^  >  jil—ik 

k= 2  i\  —ik 

gives  E[ji...L]  =  ELz(1  ~  2 l~k)x'k •  Substitution  for  x'k  from  above,  and  switching  the 
summations,  then  gives  the  required  result  in  (B5)  for  k  —  L. 

Substitution  of  (B4)  into  (B5)  and  the  result  into  the  expectation  of  (B3)  gives 
P  =  ^  ^(  — 1)4(1  —  21  *)  ^E(~  l)r  X]  Eh-ikji-jrs(\k]it) 

k=2t=2  r- o  n-4ii-ir? 

h. •••>** 

-J](-l)t(l-21-<)T(t) 

t- 2 
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where,  having  used  (b)  of  the  CL, 

n«)  =  E(-1  r  E  ci,..^»(AM,i)E(-i)‘Cl*)- 

r—t  ii---ir  k=t 

Now,  since  Efc=t(-1)*(rl*)  =  only  the  r  =  t  term  survives  in  T{t)  and  so 

T(t)  =  'y  '  Ci1...itXi1  •  ■  ■  X it 

where  the  substitution  s(X^,t)  =  A^  ■••A it  has  been  made.  This  then  reproduces  the 
expression  (9)  from  section  4.1. 


45 


DSTO-RR-0176 


46 


DISTRIBUTION  LIST 


Combinatorics  of  Crossing  Networks 
Geoff  A.  Latham 


DEFENCE  ORGANISATION 


Number  of  Copies 


S&T  Program 

Chief  Defence  Scientist  'j 

FAS  Science  Policy  > 

AS  Science  Corporate  Management  J 
Director  General  Science  Policy  Development 
Counsellor,  Defence  Science,  London 
Counsellor,  Defence  Science,  Washington 
Scientific  Adviser  to  MRDC,  Thailand 
Scientific  Adviser  Policy  and  Command 
Navy  Scientific  Adviser 
Scientific  Adviser,  Army 
Air  Force  Scientific  Adviser 
Director  Trials 

Aeronautical  and  Maritime  Research  Laboratory 

Director,  Aeronautical  and  Maritime  Research  Laboratory 

Electronics  and  Surveillance  Research  Laboratory 
Director,  Electronics  and  Surveillance  Research  Laboratory 
Chief,  Communications  Division 
Chief,  Information  Technology  Division 
Research  Leader,  Military  Information  Networks,  CD 
Head,  Network  Requirements,  CD 
Research  Leader,  Secure  Communications,  CD 
Head,  Intelligent  Networks,  CD 
Head,  Information  Access,  CD 
Head,  C3I  Networks,  CD 
Head,  Crypto-Mathematics  Research,  CD 
Head,  Information  Management  and  Fusion,  ITD 
Mr.  John  Ross,  CD  (Task  Manager) 

Dr.  Geoff  Latham,  CD  (Author) 

Dr.  Martin  Oxenham,  SSD 

DSTO  Research  Library  and  Archives 
Library  Fishermans  Bend 


1 

1 

Doc  Data  Sht 
Doc  Data  Sht 
Doc  Data  Sht 
1 

Doc  Data  Sht 
Doc  Data  Sht 
1 
1 

1 

Doc  Data  Sht 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
5 
1 

1 


Library  Maribyrnong 
Library  Salisbury 
Australian  Archives 
Library,  MOD,  Pyrmont 
US  Defense  Technical  Information  Center 
UK  Defence  Research  Information  Centre 
Canada  Defence  Scientific  Information  Service 
NZ  Defence  Information  Centre 
National  Library  of  Australia 
Capability  Systems  Staff 

Director  General  Maritime  Development 
Director  General  Land  Development 
Director  General  C3I  Development 
Director  General  Aerospace  Development 
Navy 

SO  (Science),  Director  of  Naval  Warfare,  Maritime  Headquarters 
Annex,  Garden  Island 

Army 

ASNSO  ABCA,  Puckapunyal 

SO(Science),  DJFHQ(L),  MILPO,  Enoggera,  Queensland  4057 
Air  Force 

Intelligence  Program 

DGSCD,  Defence  Signals  Directorate 
DGSTA,  Defence  Intelligence  Organisation 
Manager,  Information  Centre,  Defence  Intelligence  Organisa¬ 
tion 

Acquisitions  Program 
Corporate  Support  Program 

Officer  in  Charge,  TRS,  Defence  Regional  Library,  Canberra 

UNIVERSITIES  AND  COLLEGES 

Australian  Defence  Force  Academy  Library 

Head  of  Aerospace  and  Mechanical  Engineering,  ADFA 

Deakin  University  Library,  Serials  Section  (M  List) 

Hargrave  Library,  Monash  University 
Librarian,  Flinders  University 


1 

2 

1 

Doc  Data  Sht 
2 
2 
1 
1 
1 

Doc  Data  Sht 
Doc  Data  Sht 
Doc  Data  Sht 
Doc  Data  Sht 

Doc  Data  Sht 

4 

Doc  Data  Sht 


1 

1 

1 


1 


1 

1 

1 

Doc  Data  Sht 
1 


OTHER  ORGANISATIONS 


NASA  (Canberra) 

Info  Australia 

State  Library  of  South  Australia 
Parliamentary  Library  of  South  Australia 

CSIRO,  Mathematics  and  Information  Sciences,  North  Ryde, 
Sydney 

Investment  Technology  Group,  Australia  Inc.,  Melbourne 

ABSTRACTING  AND  INFORMATION  ORGANISATIONS 

Library,  Chemical  Abstracts  Reference  Service 
Engineering  Societies  Library,  US 

Materials  Information,  Cambridge  Scientific  Abstracts,  US 
Documents  Librarian,  The  Center  for  Research  Libraries,  US 

INFORMATION  EXCHANGE  AGREEMENT  PARTNERS 

Acquisitions  Unit,  Science  Reference  and  Information  Service, 
UK 

Library  -  Exchange  Desk,  National  Institute  of  Standards  and 
Technology,  US 

National  Aerospace  Laboratory,  Japan 
National  Aerospace  Laboratory,  Netherlands 

SPARES 

DSTO  Salisbury  Research  Library 


Total  number  of  copies: 


Page  classification:  UNCLASSIFIED 


DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
DOCUMENT  CONTROL  DATA 


1.  CAVEAT/PRIVACY  MARKING 


2.  TITLE  3.  SECURITY  CLASSIFICATION 

Combinatorics  of  Crossing  Networks  Document  (U) 

Title  (U) 

Abstract _ (U) _ 

4.  AUTHOR(S)  5-  CORPORATE  AUTHOR 

Geoff  A.  Latham  Electronics  and  Surveillance  Research  Laboratory 

PO  Box  1500 

Salisbury,  South  Australia,  Australia  5108 _ 

6a.  DSTO  NUMBER  6b.  AR  NUMBER  6c.  TYPE  OF  REPORT  7.  DOCUMENT  DATE 

DSTO-RR-  0176  AR-011-477  _  Research  Report _  June  2000 _ 

8.  FILE  NUMBER  9.  TASK  NUMBER  10.  SPONSOR  11.  No  OF  PAGES  12.  No  OF  REFS 

E8730/16/34  DST  98/224  CDS  _  45  _ |_7 _ 

13.  URL  on  the  Worldwide  Web  14.  RELEASE  AUTHORITY 

http://www.dsto.defence.gov.au/corporate/  Chief,  Communications  Division 

reports/DSTO- RR-0176.pdf _ _ _ _ 

15.  SECONDARY  RELEASE  STATEMENT  OF  THIS  DOCUMENT 
Approved  For  Public  Release 

OVERSEAS  ENQUIRIES  OUTSIDE  STATED  LIMITATIONS  SHOULD  BE  REFERRED  THROUGH  DOCUMENT  EXCHANGE,  PO  BOX  1500, 
SALISBURY,  SOUTH  AUSTRALIA  5108  _ - — — 


16.  DELIBERATE  ANNOUNCEMENT 

No  Limitations 

17.  CITATION  IN  OTHER  DOCUMENTS 

No  Limitations 

18.  DEFTEST  DESCRIPTORS 

Combinatorial  analysis 

Exchanging 

Performance  prediction 

Resources  management 

19.  ABSTRACT 

In  scenarios  such  as  software  agents  sharing  time  slots  or  blocks  in  different  physical  CPUs  or  memories, 
communications  agents  accessing  channels  or  bandwidth  on  shared  links  or  bearers  and  economic  agents 
trading  in  commodities  in  financial  markets,  the  agents  in  question  can  exchange  resources  among 
themselves  for  mutual  benefit.  This  report  develops  a  combinatorial  probability  model  to  describe  the 
like-for-like  exchange  of  resources  between  multiple  unbiased  agents.  The  expected  exchange  rates  are 
computed  for  individual  agents,  syndicates  of  agents  and  the  collection  of  all  agents.  These  performance 
benchmarks  provide  intuition  and  understanding  for  the  performance  of  real  crossing  networks.  A 
number  of  combinatorial  identities  are  also  produced  as  a  consequence  of  the  modelling  and  analysis. 


Page  classification:  UNCLASSIFIED 


